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Pre-face 


The  purpose  of  this  thesis  is  to  study  the  influence  of 
Prandtl  number  on  multicel lular  flow  in  a  narrow  vertical  slot. 
Research  such  as  this  has  practical  applications  in  technology 
systems  that  cool  or  insulate  critical  components.  Some 
examples  include  nuclear  reactor  design,  materials  processing, 
aircraft  cabin  insulation,  and  actively  cooled  aircraft  struc¬ 
tures. 

Both  low  Prandtl  number  <  Pr  <  0. 1  )  and  high  Prandtl 
number  (  Pr  >  10  )  fluids  were  investigated;  including,  Prandtl 
number  of  order  one  and  the  limiting  condition  of  Pr  0.  Due 
to  lack  of  time,  the  solution  for  Pr  -»  oo  was  not  accomplished. 
Only  the  analytical  and  numerical  equations  governing  this  type 
of  fluid  behavior  were  formulated  (  and  should  be  investigated 
in  any  follow-on  study  ). 

In  preparing  this  thesis,  I  am  eternally  grateful  for  the 
assistance  of  others  who  helped  me  through  this  project. 
Especially,  to  Capt  Daniel  B.  Fant ,  whose  superb  knowledge, 
technical  expertise,  and  patience  enabled  me  to  accomplish  my 
goal.  I  also  wish  to  thank  the  other  members  of  my  thesis 
committee,  Dr.  Franke,  and  Dr  Hitchcock.  And  finally,  thanks  to 
my  wife,  Maria,  whose  love  and  devotion  made  all  this  possible. 


David  C.  Jarman 
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Abstract 
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The  purpose  of  this  study*  is  to  investigate  the  influence 
of  Prandtl  number  on  steady,  or  unsteady,  flow  behavior  in  a 
narrow  vertical  slot.  The  focus  is  to  examine  the  effects  of 
Prandtl  number  on  natural  convective  flows,  and  the  formation  of 
single  and  multiple  cells,  over  a  range  of  Prandtl  numbers  from 

I  r  ** 1  ►. 

Pr  =  0  to  Pr  Both  numerical  as  well  as  analytical 

approaches  are  undertaken  in  studying  these  phenomena.' 

By  assuming  a  large  Rayleigh  number  and  small  gap,  the 
two-dimensional  Navier — Stokes  equations  are  reduced  to  the 
simpler  boundary-1 ayer  form.  Boundary-1 ayer  equations  are  then 

**  'IT  —  , 

derived  for  both  the  low  (  Pr''<~0.  1  )  and  high  (  Pr  St  10  ) 

Prandtl  number  regimes,  including  the  limiting  conditions  of 

>  f  1  _ 

Pr  -*  0  and  Pr  -*\jc*.  Furthermore,  the  equations  obtained  for  the 
low  Prandtl  number  regime  also  hold  for  any  finite  Prandtl  number 

in  the  range  0  <  Pr~'£'~i.O.  \ 

y 

Numerical  solutions  to  these  equations  are  obtained  by 
using  a  stable,  second-order ,  f ul ly-impl ici t ,  time-accurate, 
Gauss-Siedel  iterative  procedure.  In  addition  to  numerical 
solutions,  steady-state  analytical  solutions  are  derived  using 
perturbative  methods.  {-  I _ 

For  most  of  the  Prandtl  numbers  studied,  the  boundary-1 ayer 
approximations  produced  results  that  compared  well  with  those 
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obtained  by  other  researchers  using  the  complete  Navier — Stokes 
equations  and  more  sophisticated  numerical  and  analytical 
techniques.  The  exception  is  in  the  high  Prandtl  number  regime, 
which  may  have  been  due  to  neglecting  the  streamwise  viscous 
terms.  However,  the  results  show  that  the  boundary— 1 ayer 
approach  captures  all  o-f  the  important  physics  of  low  Prandtl 
number  fluid  flows,  flows  where  Pr  ~  1.0,  and  a  significant 
portion  of  that  involved  with  fluids  of  high  Prandtl  number. 

Due  to  lack  of  time,  the  limiting  condition  where  Pr  -»  oo 
was  not  accomplished.  However,  the  equations  for  this  type  of 
fluid  are  developed  and  presented  in  Chapter  III.  The  Pr  -♦  oo 
case  should  be  investigated  in  any  follow-on  studies  that  are 
conducted. 


MULTICELLULAR  NATURAL  CONVECTION  IN  A  VERTICAL  SLOT 


I_.  Introduction 

The  study  of  natural  convection  in  vertical  slots  has 
received  considerable  attention  from  researchers  over  the  past 
few  decades.  Practical  applications  of  this  research  include 
nuclear  reactor  design,  aircraft  cabin  insulation,  cooling  of 
electronic  equipment,  thermal  storage  systems,  and  material 
processing.  At  present,  particular  emphasis  is  being  placed  on 
the  design  of  hypersonic  vehicle  airframe  structures  where  the 
outer  skin  is  cooled  by  the  aircraft's  fuel.  Knowledge  of  the 
boundary  conditions  that  generate  multicellular  flows  is 
critical  to  the  design  of  such  vehicles  since  the  formation  of 
convective  cells  is  the  primary  mechanism  for  heat  transfer. 

A  complete  understanding  and  analysis  of  the  flow  in  the 
boundary-1 ayer  regime,  corresponding  to  the  limiting  conditions 
of  high  Rayleigh  number  and  small  gap  width,  is  still  lacking. 
In  addition,  the  influence  of  Prandtl  number  on  the  flow  field 
is  not  well  understood.  In  most  cases,  the  focus  of  previous 
investigations  has  been  on  numerical  studies  of  fluids  with  low 
Rayleigh  numbers  or  high  Prandtl  numbers  <  Pr  >  10.0  ). 

This  thesis  investigates  the  steady,  or  unsteady,  flow 
behavior  in  a  narrow  vertical  slot  induced  solely  by  buoyancy 
effects.  The  intent  of  this  study  is  to  examine  the  character- 
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istics  of  natural  convective  -flows  over  a  range  of  Prandtl 
numbers  from  Pr  =  0  to  Pr  =  oo  —  including  low  Prandtl  numbers 
(  Pr  <  0.  1  ),  Prandtl  number  of  order  one  (  Pr  ~  1.0  ),  and  high 
Prandtl  numbers  (  Pr  >  10.0  )  —  and  how  the  nature  of  the 
flow  instability  changes  with  Prandtl  number. 

In  low  Prandtl  number  fluids  thermal  diffusivity  dominates, 
therefore  the  instabilites  are  hydrodynamic  in  nature.  Whereas, 
in  high  Prandtl  number  fluids  momentum  diffusivity  dominates  and 
the  instabilities  are  more  thermal  in  nature. 

Two  approaches  are  used  to  accomplish  the  objectives  of 
this  thesis.  First,  the  2-D  Navier — Stokes  equations  are 
simplified  for  the  conditions  of  large  Rayleigh  number  and  small 
gap  width.  This  reduces  the  equations  to  a  2-D  boundary-1 ayer 
form.  Finite-differencing  techniques  are  used  to  discretize  the 
equations,  which  are  then  solved  by  computer  (  Vax  8650  )  to 
obtain  the  stream  function,  vorticity,  and  temperature 
distributions  within  the  vertical  slot. 

It  should  be  noted,  there  was  considerable  doubt  that 
the  boundary-layer  form  of  the  Navier — Stokes  equations  would 
capture  the  instabilities  present  in  natural  convective  flow. 
This  is  because  the  streamwise  viscous  terms  are  neglected  and 
it  was  believed  that  these  terms  play  an  important  role  in 
generating  the  multicellular  instability.  Fant  (  1987  ) 
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demonstrated  that  this  is  not  true  -for  flow  through  an  annulus 


with  a  small  radial  gap.  All  of  the  physical  properties  of  the 
flow  were  reproduced  and  compared  well  with  experimental  data. 
The  same  should  also  be  true  for  the  vertical  slot  as  long  as 
the  gap  is  narrow  and  the  Rayleigh  number  is  large. 

Second,  perturbation  methods  are  used  to  determine  the 
analytical  steady-state  solution  to  the  boundary-layer 
equations.  These  solutions  are  then  compared  to  the  numerical 
results  and  serve  to  reinforce  confidence  in  the  numerical 
algorithm  and  procedure. 

Li terature  Review 

Before  discussing  the  details  of  this  thesis,  a  review  of 
related  literature  in  this  field  of  study  is  appropriate. 

De  Vahl  Davis  and  Mallinson  (  1975  )  studied  natural 
convection  in  a  vertical  slot  for  a  fluid  with  Pr  =  lOOO  by 
numerically  solving  the  conservative,  second-order ,  finite- 
difference  approximations  to  the  steady-state  Boussinesq 
equations.  The  left  vertical  wall  was  hot,  while  the  right  one 
was  cool.  The  flow  was  assumed  to  be  skew-symetric  about  the 
horizontal  centerline  and  no  attempt  was  made  to  perturb  the 
flow.  Random  occurrence  of  round-off  errors  in  the  computations 
was  sufficient  to  initiate  secondary  motion.  Several  stream 
function  and  temperature  plots  were  generated  for  increasing 
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values  of  Rayleigh  number.  Their  results  show  that  a  tempera¬ 
ture  inversion  within  the  -flow  -field  occurs  when  convective 
cells  are  formed  in  fluids  with  high  Prandtl  numbers. 

Lee  and  Korpela  (  1983  )  studied  the  behavior  of  fluids  in 
a  narrow  vertical  slot.  Their  success  in  numerically  resolving 
the  multicellular  flow  field  was  attributed  to  the  fourth-order 
differencing  of  Arakawa  (  1966  ),  which  was  used  to  approximate 
the  nonlinear,  convective  terms.  However,  they  solved  the 
complete  Navier — Stokes  system  of  equations: 

S-S*  «<.*>  -£♦  ^ c 

M  -  g  J  <  T,*  >  -  ^  T2  T 
T2  *  *  -  C 

where  6  is  the  Grashof  number,  Pr  is  the  Prandtl  number,  H  is 
the  aspect  ratio,  T  is  the  temperature,  J  is  the  Jacobian,  C  is 
the  vorticity,  and  is  the  stream  function.  Steady,  multi¬ 
cellular  flow  for  Prandtl  numbers  ranging  from  0  to  1000  was 
obtained.  However,  as  Prandtl  number  increased,  the  aspect 
ratio  had  to  be  increased  to  trigger  the  formation  of  multiple 
cells. 


Chait  and  Korpela  <  1989  )  numerically  simulated  the 
multicellular  flow  between  two  vertical  parallel  plates  using  a 
time-spl i ting ,  pseud ospectral  method.  The  spectral  approx i- 
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mation  was  obtained  by  expressing  the  -flow  variables  as  a  double 
sum: 

K  N 

Q  (x,y,t  )  =  Re  £  ^  ^  Q  (  n,k,t  )  (x)  ^niky/X  j  ^ 

k  =  — K+l  n=0 

in  which  Re  denotes  the  real  part  of  the  expression  that 
follows.  This  representation  uses  a  spectral  sum  of  Chebyshev 
polynomials,  T  (X) ,  in  X  and  a  Fourier  expansion  in  Y.  The 
parameter  X  is  the  non-dimensional  wavelength  of  the  multi¬ 
cellular  structure.  The  steady  flow  of  air  (  Pr  =  0.7  )  and 
time-periodic  flow  of  oil  <  Pr  =  1000  )  were  investigated. 

The  authors  stated  that  the  spectral  transfer  of  kinetic  energy 
and  thermal  variance  for  air  was  smooth  and  dominated  by  viscous 
forces;  while  that  for  oil  was  bumpier. 

Henkes  and  Hoogendoorn  (  1989  )  submitted  a  paper  this  year 
on  a  numerical  study  of  laminar  natural  convection  in  a  square 
cavity  for  air  (  Pr  =  0.71  )  and  water  (  Pr  =  7.0  )  using  the 
unsteady,  2-D  Navier-Stokes  equations.  The  cavity  consisted  of 
a  hot  vertical  wall  and  either  conducting  or  adiabatic  hori¬ 
zontal  walls.  In  the  case  of  conducting  walls,  stability  was 
lost  at  a  lower  Rayleigh  number  than  that  for  the  case  of 
adiabatic  walls.  The  flow  of  water  was  found  to  be  more  stable 
than  that  of  air.  Above  the  critical  Rayleigh  number,  the 
unsteady  flow  had  a  single  frequency  oscillation;  except  for 
air,  which  displayed  two  distinct  frequencies  when  both  the 
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horizontal  walls  were  adiabatic.  The  authors  stated  that  the 
instability  seemed  to  be  related  to  a  hydraul i c- jump  instabil¬ 
ity.  Whereas,  for  the  case  of  conducting  and  adiabatic 
horizontal  walls,  the  instablility  seemed  to  be  related  to  a 
Rayleigh-Taylor/B6nard  and  Tol lmi en-Schl ichting  instability. 

Keyhani ,  Prasad,  and  Cox  (  1988  )  performed  an  experimental 
study  of  natural  convection  heat  transfer  in  a  vertical  cavity 
with  discrete  heat  sources.  One  vertical  wall  was  isothermal ly 
cooled.  The  opposing  wall  had  eleven  alternately  unheated  and 
flush-heated  sections  of  equal  height.  The  cavity  was  filled 
with  ethylene  glycol  <  Pr  ~  20  )  and  flow  visualization  pictures 
were  obtained  showing  the  primary,  secondary,  and  tertiary  flow 
patterns.  Their  study  indicated  that  the  local  Nusselt  numbers 
obtained  in  the  enclosure  are  nearly  the  same  as  those  reported 
for  a  wide  flush-mounted  heater  on  a  vertical  plate. 

Elder  (  1965  )  performed  an  experimental  study  of  natural 
convection  for  a  liquid  (  Pr  =  1000  )  in  vertical  slots  with 
isothermal  walls.  Single  cell  flows  resulted  when  the  Rayleigh 
number  was  less  than  10^.  For  larger  aspect  ratios,  a  steady 
secondary  flow  was  superimposed  on  the  basic  single  cell. 
Increasing  the  Rayleigh  number  to  above  10^*,  produced  tertiary 
flow  with  counter — rotation. 

Elder  (  1966  )  numerically  solved  the  same  problem 
described  above  for  a  moderate-size  vertical  slot.  He  dupli- 
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cated  the  basic  flow  field  obtained  in  the  experiment,  bat  was 
unable  to  capture  the  secondary  flow.  Elder  stated  that  the 
nonlinear  terms  began  to  dominate  the  motion  as  the  Rayleigh 
number  was  increased;  and  was  able  to  numerically  show  the 
development  of  the  boundary-1 ayers  and  fully  developed 
boundary— 1 ayer  flow. 

Vest  and  Arpaci  (  1969  )  performed  an  analytical  investi¬ 
gation  of  natural  convection  in  a  narrow  vertical  slot  using 
linearized  hydrodynamic  stability  theory.  They  obtained  a 
neutral  stability  curve  for  the  conduction-dominated  flow 
regime.  In  addition,  stream  function  plots  for  the  secondary 
flow  fields  were  obtained  that  qualitatively  agreed  with 
experimental  results. 

Korpela  et  al^.  <  1973  )  used  linear  stability  theory  to 
examine  the  stability  of  the  conduction  regime  for  natural 
convection  in  a  vertical  slot.  They  claimed  that  instability 
set  in  as  horizontal  cells  for  Pr  <  12.7  .  The  instability  was 
thought  to  be  hydrodynamic  in  origin,  resulting  from  the 
vorticity  distribution  of  the  base  flow. 

Korpela  (  1974  )  studied  the  problem  of  the  narrow  slot 
inclined  at  an  angle  6,  from  the  vertical.  Instabilities  set  in 
as  transverse,  traveling  waves  for  small  angles  of  inclination; 
whereas  longitudinal  cells  formed  as  6  reached  a  certain  value. 
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For  Pr  <  0.24,  Korpela  claimed  that  only  horizontal  cells  were 
possible  and  that  the  stability  of  the  flow  was  mainly  a 
function  of  Pr  tan  6. 

Pepper  and  Harris  <  1977  )  numerically  obtained  2-D  natural 
convective  flow  patterns  in  rectangular  and  annular  vertical 
cavities.  The  energy  and  vorticity  equations,  in  divergence 
form,  were  solved  using  central-dif f erencing  and  a  strongly 
implicit  procedure.  A  weak  multicellur  flow  pattern  was 
obtained  for  Pr  =  1000  and  Ra  =  5  x  105  . 

Seki  et  al^.  (  1978  )  performed  an  experimental  analysis  of 
natural  convection  in  narrow  vertical  rectangular  cavities. 
Transformer  oil,  water,  and  glycerin  were  used  as  working 
fluids.  Multicellular  secondary  flow  was  observed  in  each  case. 
In  the  case  for  oil,  as  the  temperature  was  increased,  tertiary 
flow  with  counter — rotating  cells  developed  until  the  flow  near 
the  upper  region  of  the  hot  wall  became  unsteady  and  turbulent. 
They  concluded  that  the  flow  more  easily  shifted  from  laminar  to 
turbulent  as  the  Prandtl  number  decreased  and  the  cavity  width 
i ncreased. 

From  the  literature  search,  it  is  evident  that  there  is 
considerable  interest  in  studying  the  multicellular  flow 
phenomenon.  Many  theoretical,  numerical,  and  experimental 
studies  have  been  accomplished.  However,  none  have  used  the 
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boundary-layer  approach  to  solving  the  Navier — Stokes  equations 
for  the  vertical  slot  geometry. 
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This  chapter  describes  a  vortici ty-stream  function 
formulation  of  the  Navier — Stokes  equations  in  dimensionless 
form.  Using  this  approach,  the  pressure  terms  are  eliminated 
and  the  conservation  of  mass  principle  is  automatically 
satisfied.  The  solutions  to  these  equations  are  used  to 
generate  stream  function,  vorticity,  and  temperature  contours 
that  are  useful  for  visualizing  and  analyzing  the  flow  field. 
The  final  form  of  the  equations  is  presented  in  Eqs  (  2.30  - 
2.32  )  and  (  2.44  -  2.46  ). 

Physical  Model 

a.  The  flow  is  unsteady  and  two-dimensional 
(see  Figure  1 ) . 

b.  The  fluid  is  initially  at  rest. 

c.  The  left  vertical  wall  is  isothermally  heated,  the 
right  vertical  wall  is  isothermally  cooled. 

d.  The  fluid  is  Newtonian  and  laminar  fluid  motion  is 
induced  solely  by  buoyancy  effects. 

e.  All  material  properties  are  constant  and  density 
variations  are  allowed  to  occur  by  means  of  the 
Boussinesq  approximation. 
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Dimensional  Formulation 


Governing  Equations  Usi nq  the  Vorti ci ty-Str earn  Function 
Approach .  The  vorti ci ty-transport  equation  is  obtained  by  taking 
the  curl  of  the  Navi er-Stokes  equations  (  ref  Incompressible 
FI ow  ,  p  328  ).  This  equation  can  be  written  as 

SU  —  —  —  _  _  —  —■?  —  _  _ 

—  +  <V-7)  W  =  <W-7>  V  +  v  v  W  +  7  x  F  ,  (  2.1  ) 

at  B 

where  the  bar  over  the  variable  indicates  a  dimensional 
quantity.  W  and  V  represent  the  dimensional  -form  of  the 
vorti city  and  velocity  vectors,  respectively. 

The  left  side  of  the  equation  is  the  total  rate  of  change 

of  particle  vorticity.  The  terr  <WV>  V  is  identically  zero  for 

2— D  flows  because  the  vorticity  vector  is  always  perpendi cul ar 

to  the  plane  of  flow.  The  term  v  T2  W  is  the  net  rate  of 

vorticity  diffusion  due  to  viscosity.  The  last  term,  7  x  F  ,  is 

o 

the  rate  of  inertial  vorticity  generated  by  the  body  forces 
which  result  from  density  perturbations  in  natural  convective 
f 1 ows . 


It  is  interesting  to  note  that  the  pressure  (  F  )  terms  do 
not  appear  explicity  in  Eq  <  2.1  ),  since  7  x  TP  =  0. 

Therefore,  the  flow  field  can  be  determined  without  knowledge  of 
the  pressure  distribution.  Using  the  Boussinesq  approximation, 
the  net  body  force  can  be  obtained.  The  resulting  expression 
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B 


=  -g/?(T  -  T  )j 


(  2.2  ) 


where  T  is  a  known  initial  temperature,  ft  is  the  coefficient  of 

o 

thermal  expansion,  and  g  is  the  acceleration  due  to  gravity. 
Making  this  substitution,  the  final  form  of  the  vorticity- 
transport  equation  becomes 


»±  +  HiigL  =  , 

at  0<x, y>  ax 

and  the  thermal  energy  and  stream  function  equations  are 


(  2.3  ) 


AT 

at 


and 


W 


=  a  T2  T 

a(x, y> 


=  T2  f 


(  2.4  ) 


(  2.5  ) 


In  the  above  equations,  X  and  Y  are  the  coordinate 
positions  in  the  f  and  j*  directions,  T  is  the  dimensional 
temperature  variable,  and  t  is  dimensional  time.  v  and  a  are 
the  momentum  and  thermal  dif fusi vities,  respecti vel y.  The 
spacial  velocities  u  and  v  have  been  written  in  terms  of  the 
stream  function,  f,  defined  as 


at 

u  =  -  — —  and 

ay 


v  = 


at_ 

ax 


(  2.6  ) 


This  reduces  the  number  of  dependent  variables  by  one.  The 
convective  terms  are  represented  by  the  Jacobian,  <><P,Q> /£(X, Y>  , 


II 


4 


where 

d(P,Q)  dP  dQ  dP  dQ 

a<x,Y>  ax  dY  aY  ax  ‘  (  2 

These  three  coupled  equations  (  Eqs  2.3  -  2.5  )  describe  the 

vorticity,  temperature,  and  stream  function  for  the  vertical 
slot. 


Boundary  Conditions.  For  temperature,  the  left  and  right 
vertical  walls  are  isothermal ly  heated  and  cooled,  respectively. 
The  top  and  bottom  horizontal  walls  are  adiabatic.  Therefore, 
the  boundary  conditions  ares 


T (0, Y)  =  T  (  2.8  ) 

T  <L,Y>  =  Te  (  2.9  ) 

—  ( X ,0)  =  (X, 1 )  =  0  .  (  2.10  ) 

aY  aY 

The  velocity  components,  u  and  v,  are  identically  zero 
along  the  vertical  and  horizontal  walls.  This  requirement  gives 
the  boundary  conditions  for  stream  function  and  vorticity: 


f  (O,  Y) 


f  <L, Y)  = 

W 

W 


f  <X,0) 


a2? 

ax2 


a2f 

a72 


=  f(X,l) 


X=0,L 


Y=0 , 1 


0 


(  2.11  ) 

(  2. 12  ) 

(  2. 13  ) 
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Initial  Conditions.  Since  the  -fluid  is  initially  at  rest, 

7  =  W  =  0  and  T  =  T  (2.14) 

o 

throughout  the  vertical  slot,  except  at  the  left  wall  where 

T  =  T.  (2.15) 

is  enforced. 

Pi mensi onl ess  Formulation 

Writing  equations  in  dimensionless  form  has  three  disctinct 
advantages: 

a.  it  simplifies  the  mathematical  or  numerical 
expressions,  allowing  for  better  interpretation 

b.  measurement  scales  are  no  longer  an  intrinsic  part  of 
the  physical  quantities,  therefore,  any  laws  governing  physical 
variables  are  valid  for  all  different  measurement  systems 

(  Fant,  1987  ) 

c.  fewer  variables  are  used  and  the  proper  dimensionless 
groups  characterizing  a  particular  flow  model  become  evident. 

Coordinate  Transf or mat ion.  The  X  and  Y  spacial  coordinates 
are  non-dimensional i zed  using  the  appropriate  length  scales  from 
the  vertical  slot  geometry, 

X  =  -  and  Y  =  -  .  (  2.  16  ) 

L  1 


II 


6 


Time,  the  other  independent  variable,  is  non-di mensi onal i zed  as 


t  = 


tv 

~2 


or 


1*  1 
The  dimensionless  gap,  G,  is  de-fined  as 


Fa 
2  * 


G  =  -  . 

1 


A1  so. 


and 


Ra  = 


Pr  =  - 

a 


gtfl  (T.  -  T  ) 
_ _ 1 _ ° _ 

vex 


Navier — Stokes  Equations  -for  Low  Prandtl  Number 
(  Pr  <  0.  1  )  .  In  dimensionless  form,  the  dependent 
are  written  as: 


T  = 


T  -  T 


T.  -  T 

l  < 


«12 


W  =  — 


(temperature) 


(vortici ty) 


In  addition. 


f  -  - 


(stream  function) 


6 

ax 

a 

l  a 

ax 

ax 

ax 

l  ax 

£_ 

<JY 

a_ 

i  a 

<9Y 

av 

aw 

\aw 

(  2.17  ) 

(  2.18  ) 

(  2.19  ) 

(  2.20  ) 

variables 

(  2.21  ) 

(  2.22  ) 

(  2.23  ) 

(  2.24  ) 

(  2.25  ) 
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(  2.26  ) 


Zi 

ax2 


a2 


2  2 

l  ax^ 


a 

av2 


i  £_ 
2  2 
i  aY^ 


u 


af 

ay 


yu 

1 


v  = 


at 

ax 


yv 

L 


When  Eqs  (  2. 16  -  2.29  )  are  substituted  into  Eqs 
2.5  ),  the  -following  non-dimensional  -form  o-f  the  Navier 
equations  result 
Energy: 

ay  > 


Pr 


91  f  at  91 

at  \  9X  9V 

af  aT  1 

.  i  f 

921 

ay  ax  J 

Pr  ( 

ax2 

Vortici ty: 

f  G2  <*»  +  G  r 

aw  aw 

af  I  "l 

\  6  at  +  6  L  ax 

aY  ~  ax 

aY  J  j 

t  ax^ 

+  g2 

a2w 

aY2 

Stream  Function: 

W  = 

-cr 

♦  ^5 

— „ 

b  ax 

9V2 

91 

ax 


(  2.27  ) 

(  2.28  ) 

(  2.29  ) 

(  2.3  - 
Stokes 


(  2.30  ) 


(  2.31  ) 

(  2.32  ) 
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where 


7 


2 


(  2.33  ) 


G,  Pr ,  and  Ra  are  three  dimensionless  parameters  that  arise 
from  this  analysis  and  are  used  to  simulate  various  flow 
conditions  and  geometries.  The  gap  number,  G,  defines  the 
vertical  slot  geometry.  It  is  analagous  to  the  aspect  ratio 
used  to  define  other  geometries,  such  as  airfoils.  The  Prandtl 
number,  Pr ,  character i zes  the  fluid  and  is  the  ratio  of  momentum 
to  thermal  diffusivity.  The  Rayleigh  number,  Ra,  is  the  ratio 
of  the  buoyant  forces,  which  are  destabilizing,  to  the  viscous 
forces,  which  tend  to  stabilize  the  flow.  The  magnitude  of  Ra 
can  be  used  to  predict  the  onset  of  thermal  and/or  hydrodynamic 
instability,  since  the  transition  to  instability  usually  occurs 
at  some  critical  value  of  Ra. 


The  variables  W,  f,  and  t,  in  Eqs  (  2.30-2.32  ),  are 
non-dimensional i zed  with  u,  instead  of  ot,  to  reflect  a 
thermal-dif fusion  dominated  energy  equation  for  small  Prandtl 
numbers.  Later  on  in  this  thesis,  when  high  Prandtl  number 
flows  are  investigated,  a  will  be  used  to  non-dimensional i ze  the 
equations  to  reflect  momentum-diffusion  dominated  flows. 

Non-di mensi onal  Boundary  Conditions.  The  non-dimensional- 
ized  boundary  conditions  for  the  vertical  slot  are  obtained  by 
subtituting  Eqs  (  2.21  -  2.23  )  into  Eqs  <  2.11  -  2.13  )  and 
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• 

noting  that: 

at  X  = 

0  , 

X  = 

0  and  T  =  1  ; 

at  X  = 

L  , 

X  = 

1  and  T  =  0  ; 

at  Y  * 

0  , 

Y  = 

~  «*T 

0  and  a7  =  0  5 

at  Y  = 

1, 

Y  = 

1  and  ay  ~  °  ' 

The  boundary 

conditions 

are: 

• 

■f  <0 ,  Y) 

■  f  <1 

,Y) 

=  f (X ,0)  =  f(X,l)  =  0 

( 

2.34 

) 

d* 

ax 

x=o , 

frf 

d\  0 

1  Y=0 , 1 

( 

2.35 

) 

T  <0 , 

Y) 

=  1,  T ( 1 , Y)  =  0 

( 

2.36 

) 

O 

II 

o 

II 

> 

*1* 

( 

2.37 

) 

w  .  i_ 

( 

2.38 

> 

6  dX  X=0 , 1 

and 

■ 

II 

3 

( 

2.39 

) 

• 

9V  ' Y=0 , 1 

II  -  10 

Non-di mensi onal  Initial  Conditions.  The  dimensionless 


initial  conditions  are  similar  to  Eqs  (  2.  14  -  2.15  )  -for  the 
entire  vertical  slot, 

f  =  W  =  T  =  0  .  (2.40) 

Except,  at  the  left  vertical  wall  where 

T  =  1  .  (  2.41  ) 

Note  that  at  rest,  Ra  =  0,  and  the  problem  becomes  one  of 

steady-state  conduction.  Initial  conditions  other  than  zero  are 

also  used,  and  will  be  explained  in  Chapter  IV. 

Navi er -Stokes  Equations  -for  High  Prandtl  Number 
(  Pr  >  10.0  ) .  The  procedure  for  formulating  the  high  Prandtl 
number  equations  is  similar  to  that  for  small  Prandtl  number.  a 
replaces  i>  as  the  non-dimensional izing  parameter  to  reflect  a 
momentum— di ffusi on  dominated  vorticity  equation.  Dimensionless 
vorticity,  stream  function,  and  time  are 

-  2  -  - 
■  i  W  1  ,  f  ,  .  t  a  .  _  . 

W  =  -  ,  f  =  —  ,  and  t  =  — =■  .  (  2.43  ) 

a  *  a  1  j2 

The  dimensionless  form  of  the  energy,  vorticity,  and  stream 

function  equations  become 
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Energy: 


G 


2  3T 

at 


G 


f  £f  aT 
{  ax  aw 


af  ar  \ 
aw  ax  J 


a 2t  2  a2 t 

—  +  G  — 2 

ax  ay 


(  2.44  ) 


Vorti ci ty: 


„2  aw 
G  at 


at  aw 
ax  ay 


<m  aw  "I 
ay  ax  j 


=  Pr 


{ 


__  ar 
Rai  ax 


a2  w 
ax2 


G 


2  a"w  1 


(  2.45  ) 


Stream  Function: 


W 


2  2 

i  a_f  at 

2  2  2 

g  ax  ay 


(  2.46  ) 


Compare  these  equations  with  Eqs  (  2.30  -  2.32  ).  The 
energy  equation  is  no  longer  dependent  on  Prandtl  number;  but  a 
strong  dependency  now  exists  on  the  diffusion  side  of  the 
vorti city  equation.  Boundary  and  initial  conditions  are  the 
same  (  Eqs  2.34  -  2.41  ). 
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Ill 


Analytical  Analysi s 

In  this  chapter,  simplified  boundary-1 ayer  equations  are 
derived  from  the  2-D  Navi er-Stokes  equations  for  certain 
limiting  conditions.  Analytical  steady-state  solutions  to  these 
equations  are  then  determined.  The  best  approach  to  this 
analysis  was  taken  from  Fant  (  1987  ). 

Boundary-Layer  Equations  for  Low  Prandtl  Number  (  Pr  <  0.1  )  . 

Using  a  high  Rayleigh  number /smal 1 -gap  asymptotic 
expansion,  Fant  (  1987  )  derived  simplified,  boundary-1 ayer 
equations  from  the  2-D  Navier — Stokes  equations  that  modeled  a 
narrow  cylindri  d  annulus.  These  equations  were  valid  in  the 
double  limit  £.=  Ra  -♦  oo  and  G  -*  0  .  The  following  expansions  for 
the  dimensionless  velocity  components,  vorticity,  stream 


function,  and  temperature  were  obtained: 

1/4  ~  -1/4 

u  =  Ra  u  +  0 (Ra  >  (  3.1  ) 

v  =  Ra1/4  v  +  0(1)  (  3.2  ) 

W  =  Ra374  5J  +  0(Ra1/4)  (  3.3  ) 

f  =  Ra1/4  ?  +  0(Ra~1/4)  (  3.4  ) 

T  =  T  +  0(Ra”1/2>  .  (  3.5  ) 


The  tilda  over  the  variable  is  used  to  indicate  a  scaled 
dimensionless  quantity.  These  scalings  are  a  result  of  the 
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boundary  layer  thickness  being  of  order 

6  ~  Ra_1/4  .  (  3.6  ) 

The  above  scalings  will  be  incorporated  into  the  present 
study  to  obtain  boundary-1 ayer  equations  that  model  the  narrow 
vertical  slot  geometry.  Since  the  left  and  right  wall  boundary 
layers  merge  as  G  -»  0,  the  point  at  which  they  just  begin  to 
contact  each  other  is 

26  =  L  .  (  3.7  ) 

Substituting  Eq  (  2.18  )  into  Eq  (  3.7  )  gives  the 
relationship  between  gap  number  and  boundary-1 ayer  thickness. 

26  =  G1  or  G  =  -  6  .  <  3.8  > 

1 

-1/4  -1/4 

But,  6  ~  Ra  .  Therefore,  G  ~  Ra  .  With  this  result, 

the  proper  scale  for  gap  number  is 

~  -1/4 

G  =  G  Ra.  (  3.9  ) 

To  maintain  the  unsteady  terms  in  Eqs  (  2.30  -  2.31  ),  the 
appropriate  scale  for  time  is 

~  -1/2 

t  =  t  Ra  .  (3.10) 

Eqs  (  3.3  -  3.5  )  and  (  3.9  -3.10  )  can  now  be  subtituted 
into  Eqs  (  2.30  -  2.32  ).  Starting  with  the  energy  equation 
(  Eq  2.30  ) , 
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<Ra1/2>  *1 


_  f  Ra  1/4  a?  £T  _  1/4  a?  dT  \ 

1/4  |  ax  aw  aw  ax  j 


Pr  l  Ra1/2  aY2  /  ’ 


which  reduces  to 


G2  *1 


^  ~  ~  ^  <V  v 

f  di  dJ  _  di  6T  ] 

\  ax  aw  aw  ax  j 


-  i  f  *2t  ,  g2  a2T  | 

Pr  l  ax2  Ra1/2  aY2  /  ' 


(  3.11  ) 


Likewise,  substituting  into  the  vorticity  equation  (  Eq  2.31  ) 


gives 


Pr  G2  **  +  g  pr  /  —  —  -  ™  \ 

~  \  ax  aw  ax  aw  J 


G 


f  a2 w  g2  a2w  \ 

l  ax2  Ra1/2  aw2  J 


(  3.12  ) 


And,  the  stream  function  equation  (  Eq  2.32  )  becomes 


Ra1/2  av2 


~2  'v 

G  W 


(  3.13  ) 


For  large  Rayleigh  number  and  small  gap,  the  following 
orders  of  magnitude  are  evident: 


eP _  a2? 

Ra1/2  aw2 


Ill 


3 


2~ 

aw 


Ra1/2  aY2 


<< 


a2w 

ax2 


G2 


_ 

Ra1/2  aY2 


<< 


a2? 

ax2 


Neglecting  these  higher  order  terms,  the  Navi er-Stokes 
equations  -For  low  Pr andtl  numbers  reduce  to  the  simplified 
boundary-1  ayer  -form 

Energy: 


G2  *1 

at 


v  ~  'V  >Vi  - 

fatal  -  _  1 

\  ax  ay  aY  ax  )  Pr 


a^r 

ax2 


(  3.14  ) 


Vortici ty: 


Pr  G2  fl 

*** 

at 


+  G  Pr 


f  af^  aw  _  aw  at  1 
|  ax  aY  ax  aY  ] 


F  aT 

G  ax 


Pr 


a2w 

ax2 


(  3.15  ) 


Stream  Function: 


2~ 
1_  a^-F 

~2  2 

g  ax^ 


=  w 


(  3.16  ) 


Eqs  (  3.  14  -  3.16  )  are  valid  only  -for  large  Rayleigh  and 

small  gap  numbers.  They  also  hold  for  any  finite  Prandtl  number 

less  than  or  equal  to  one.  Furthermore,  there  is  no  longer  any 

2  2 

dependency  on  d  /dY  terms.  Thus,  the  number  of  terms  are 
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reduced,  as  well  as  the  boundary  and  initial  conditions  that  go 


with  them.  The  boundary  conditions  are: 

f (0, Y)  =  f(l,Y)  =  ?(X,0)  =  f(X,l) 


=  0 


a f 

ax 


x=o,  1 


&f_ 

av 


=  o 


Y=0,  1 


T  <0,Y)  =  1,  T  (1  ,  Y)  =  0 


ai 

aw 


a  * 
ax2 


Y=0, 1 


x=o,  1 


=  o 


CNrf 

=  G  W 


The  initial  conditions  are: 

•M  ru 

f  =  W  =  T  =  0  , 

except  at  the  left  vertical  wall  where 


T  =  1 


(  3.17  ) 

(  3.18  ) 

(  3.19  ) 

(  3.20  ) 

(  3.21  ) 

(  3.22  ) 

(  3.23  ) 


Steady— State  Analytical  Solution  for  Low  Prandtl  Number  . 
To  a  leading  order ,  in  the  limit  as  G  -*  0,  Eqs  (  3.  14  -  3.16  > 
reduce  to 

Energy: 


Vorticity: 


(  3.24  ) 


(  3.25  ) 
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Stream  Function 


d-f 

— =  =  0  .  (  3.26  > 

ax 

Therefore,  the  energy  equation  (  Eq  3.14  )  simplifies  to 
the  steady-state,  1-D  conduction  equation,  which  can  then  be 
integrated  to  obtain 


T  =  l  -  x 


(  3.27  ) 


for  the  boundary  conditions  of  this  problem. 


From  the  no-slip  condition,  the  vorticity  and  velocity 
components  are  identically  zero  at  the  walls.  The  convective 
terms  can  thus  be  ignored  and  the  temperature  distribution  is 
completely  determined  by  the  linear  relationship  of  Eq  (  3.27  ). 
With  these  results,  and  assuming  steady-state  (  a/ at  ) 
conditions,  the  following  asympotic  expansions  are  obtained 
(  Fant,  1987  >: 

T  =  (  1  -  X  )  +  G4T1  +  G8T2  +  0(G16)  (  3.28  ) 

W  =  GWj  +  G5W2  +  0(G9)  (  3.29  ) 

f  =  G3^  +  G7f 2  +  0(G9)  .  (  3.30  ) 

Rearranging  the  stream  function  equation  for  low  Pr  (  Eq 
3.16  )  results  in  the  following  expression  for  vorticity: 


(  3.31  ) 
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Assuming  steady-state  conditions  and  substituting 


Eq  (  3.  31  )  into  Eq  (  3. 15  )  now  gives 

-3"  L 

ax  ax^ay  ax"  ay 


-  Pr 


f  £1  *  j  _  *_£  £1  | 

l  ax  ax2ay  ax3  aY  J 


G  51 

ax 


+  Pr 


a4f 

ax4 


(  3.32  ) 


Substituting  Eq  (  3.30  >  into  this  expression  and  equating 
like  coef -f i ci ents  yields 


a4f 


ax 


Pr 


{“ 


l 

Pr 


a3f 


ax  ax  aY 


ax3  ay 


a4f . 


ax 


} 


(  3.33  ) 


at. 


ax 


=  o 


(  3.34  ) 


~9 

6  :  Pr 


i*  ax  ay  ax 


»3*2 


ax  ax  ay 


j3<2  ^ 

ax3  ay 


G13:  Pr 


*3fl  *2 

1  * 

dJ2 

ax3  ay 

; 

ax 

*2 

*2 

ax 

s 

N 

X 

<0 

ax3 

ay 

=  o 


-  0  . 

ay  J 


<  3.35  ) 


(  3.36  ) 


An  expression  for  fj(X),  in  terms  of  the  Prandtl  number, 
can  be  obtained  by  integrating  Eq  (  3.33  ). 


Ill  -  7 


X 


x~ 


X 


VX)  =  4TpF 


+  31" 


+  C^X  +  . 


(  3.37  ) 


Applying  the  boundary  conditions  (  Eqs  3.17  -3.18  )  yields 


the  -final  result  -for  -f^(X)  . 


VX>  =  24TF 


2  2 
X  (  X  -  1  ) 


(  3.38  ) 


Equation  <  3.38  )  can  be  differentiated  twice  to  obtain  an 


expression  for  W^(X) 


a2f 


W^X)  = 


ax 


1  _  (  xx  -  X  +  —  ) 

2  Pr  1  *  *  6  1 


(  3.39  ) 


Using  Eq  (  3.34  ),  and  noting  the  df^/dV  =  0  and  dT^/dX  ~  0 
leads  to 


f2  =  w2  =  °  • 


Substituting  Eq  (  3.38  )  into  Eq  (  3.30  )  and  Eq  (  3.39  ) 
into  Eq  (  3.29  )  gives  the  final  result  for  stream  function  and 
vorticity.  ~ 


?  =  5TTF  X2  <  x  -  1  >2 


(  3.40  ) 


W  = 


2  Pr 


{ 


-  X 


H 


(  3.41  ) 


The  temperature  is  (  see  Fant,  1987  ) 


8 


T  =  (  1  -  X  )  +  G  T_  , 


(  3.42  ) 


where 


-1 


96  Pr 


h  *  >} 


34560 


(  3.43  ) 


-X 


11 


,10 


-X 


415800 


75600 


45360 


10800 


—  X 

6900 


-X 


15120 


387701 


(  3.44  ) 


,11 


-X 


10 


-X 


8 


83160 


15120 


6480 


5040 


7560 


-X‘ 


-X 


-4  YU 
55  * 


15120 


2  Y10 
5 


30240 


=5  X9 
6  * 


15 


Zl  x5 

10  * 


*  5 


997934 


3  xa 

4  * 


-X 

171 


(  3.45  ) 


-5  y7 
21 


(  3.46  ) 


It  is  important  to  note  that  the  above  analytical  expres¬ 
sions  are  valid  for  any  finite  Prandtl  number  less  than  or  equal 
to  one. 


Limiting  Equations  for  Zero  Prandtl  Number.  One  of  the 
primary  objectives  of  this  thesis  is  to  study  the  behavior  of 
very  low  Prandtl  number  fluids,  which  is  characteristic  of 
liquid  metals,  in  a  vertical  slot.  To  accomplish  this,  the  low 
Prandtl  number  boundary-1 ayer  equations  are  simplified  for  the 
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limiting  condition  of  Pr  -»  0.  Although  numerical  studies  of 
such  flows  have  been  done  before,  none  have  used  the  boundary- 
layer  equation  approach. 


First,  let  Pr  -»  0  and  assume  the  temperature  is  0(1). 

Next,  assume  that  the  vorticity  and  stream  functions  have  the 
form 

W  =  X  W  +  .  .  .  (  3.47  ) 

f  =  X2F  ■*■••••  (  3.48  ) 

And,  from  the  stream  function  equation  (  Eq  3.16  ), 

X*2  ~  B2  X t  (  3.49  ) 

By  substituting  these  expressions  into  Eqs  (  3.14  -  3.16  ), 
the  energy  equation  can  be  written  as 


G2  *1 

at 


+ 


G  X, 


f  &F  &T 

\  ax  aw 


aF  £  1 
ay  ax  j 


(  3.50  ) 


Likewise,  the  vorticity  equation  becomes 


Pr  6  Xt 


aw 


at 


Pr  G  X2 


f  aF  aw 
\  ax  ay 


aF  aw  \ 
ay  ax  j 


Pr  X„ 


(  3.51  ) 
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And,  the  stream  -function  equation  is 


W 


‘  x.  i* 
2  2 
*  ax 


(  3.52  ) 


To  retain  the  physics  of  the  flow,  it  becomes  evident  that 


G2  Pr  X. 


Pr  G  X2  X1 


Pr  Xj  ~  G  . 


(  3.53  ) 


By  equating  the  first  and  third  terms, 


~2 

G  Pr  X, 


-  Pr  Xj  , 


(  3.54  ) 


it  is  apparent  that 


t  ~  8  . 


(  3.55  ) 


Balancing  the  viscous  and  buoyancy  terms  now  gives: 

x  G 

X1  ~  Pr  ’ 

Substituting  this  expression  into  Eq  (  3.49  )  yields: 


(  3.56  ) 


^2  -  Pr  • 


(  3.57  ) 


But,  since  Pr  §  X}  X2  ~  S  <  Eq  3.51  ), 

"■[*0 (£)  ~  § 


,  §  ~  Pr' 


(  3.58  ) 


(  3.59  ) 
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From  this  and  Eq  (  3.55  ),  it  can  be  seen  that 


t  ~  Pr 


1/2 


( 


Using  these  relations  and  Eqs  (  3.9  )  and  (  3. 10  )  , 
set  of  scaled  variables  can  be  defined. 

1/4 


G  =  G  Pr 


-1/4 


t  Pr 


-1/2 


=  t 


(S) 

(S) 


Also,  from  Eqs  (  3.56  ),  (  3.57  ),  and  (  3.59  ), 

-3/4 


Xj  ~  Pr 


X2  ~  Pr 


-1/4 


And  from  Eqs  (  3.5  ),  (  3.47  ),  and  (  3.48  ), 

W 


Pr'3/4W  ♦  .  .  . 


~  -1/4 

f  ~  Pr  + 


( 


T  ~  T  + 


The  scaled  variables  now  show  a  unique  dependency  on 
Pr,  G,  and  t.  Notice  also,  that  Ra/Pr  (  in  Eqs  3.61  and 
is  defined  as  the  Grashof  number.  The  energy,  vorticity, 
stream  function  equations  can  now  be  written  as 
Energy: 


Pr 


{ 


~2  OT 
»  — 

at 


[SF  <9T 

ax  ay 


ar  ai  1  \ 
ay  ax  J  j 


( 


3.60  ) 

a  new 

3.61  ) 

3.62  ) 

3.63  ) 

3.64  ) 

3.65  ) 

3.66  ) 

3.67  ) 

Ra, 

.62  ) 
and 

3.68  ) 
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Vorticity: 


AO 

G“ 


aw 

A 

at 


+ 


G 


f  £  aw 
\  ax  ay 


£  1  p  <?T 

ax  ay  j  b  ax 


(  3.69  ) 


to 


Stream  Function: 


a2F 


W  . 


(  3.70  ) 


But,  in  the  limit  as  Pr  -*  0,  the  energy  equation  simplifies 


(  3.71  ) 


which  is  now  completely  uncoupled  from  the  vorticity  equation. 
Eq  (  3.71  )  can  be  integrated  twice  to  obtain  an  expression  for 
the  temperature  distribution  as  a  function  of  X  . 


T  =  CjX 


(  3.72  ) 


By  applying  the  boundary  conditions,  T.  and  T  , 

1  ° 

T  ( X )  =  1  -  X  . 


Thus, 


and 


£T 

ax 


-l 


aT 

a?  = 


o 


(  3.73  ) 

(  3.74  ) 

(  3.75  ) 


This  last  expression  automatically  satisfies  the  adiabatic 
boundary  condition  on  the  top  and  bottom  horizontal  walls. 

Using  the  above  results,  the  vorticity  equation  reduces  to  its 
final  form. 


Ill 


13 


~2  £W  +  g  f  —  \ 

~  \  ax  d'i  ax  ay  j 


=  -  G  + 


d2\* 

ax2 


The  boundary  conditions  -for  the  equations  are: 


F  (0 ,  Y)  =  F  ( 1 ,  Y)  =  F  ( X  ,  0)  =  F  <  X  ,  1 )  =  0 


dF 

ax 


X=0, 1 


£F 

0Y 


=  0 


Y=0 , 1 


T  (0,  Y)  =  1,  T  ( 1 ,  Y)  =  0 


ai 

a\ 


*2F 


=  0 


Y=0,1 


ax  1 x=o , l 


e  w  , 


and  the  initial  conditions  remain 

F  =  W  =  T  =  0  , 

except  at  the  left  vertical  wall  where 


T  =  1 


(  3.7 b  ) 

(  3.77  ) 

(  3.78  ) 

(  3.79  ) 

(  3.80  ) 

(  3.81  ) 

(  3.82  ) 

(  3.83  ) 


Steady-State  Analytical  Solution  for  Pr  ■*  0  .  Rearranging 
the  stream  function  equation  (  Eq  3.70  )  for  the  limiting 
condition  of  Pr  -»  0  results  in  the  following  expression  for 
vorti ci ty : 


W 


1  <>2f 


<  3.84  ) 
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Assuming  steady-state  conditions  and  substituting  Eg 


(  3.84  )  into  Eg  (  3.7 6  )  yields 


(  d3|r  _  df_ 

l  ax  a2x  aw  ax3  aY 


(  3.85  ) 


In  the  limit  as  G  -*  0,  the  perturbative  solution  to  this 
eguation  can  be  obtained.  The  expansion  -for  the  stream  function 
i  s 

F  =  G  +  G  F^  +  0<G*  )  .  (  3.86  ) 

And,  vorticity  can  be  expressed  as 

W  =  G  Wj  +  G5  W2  +  0 (G9)  .  (  3.87  ) 

Note  that  the  effects  of  the  convective  terms  first  appear  in  F2 
and  W?.  Substituting  Eg  (3.86  )  into  Eg  (  3.85  )  gives 


{ 


7  ^1  *^1 


ill  ^2  *^1 


+  G 


11  ^1  *^2 


ax  ax*"aY 


215  ^2  **f2 
G  — - — 


ax  ax  ay 


ax  ax  aY 


27  ^1  ^1  211  *3p2  ^1 

G  ^ —  -  —  G 


-  G 


ax  axaY 


11  *3fi  ^2 

ax3  aY 


ax3  aY 


ax3  aY 


-  G 


15  **F2  * 

ax2  aY 


— } 
aY  J 


=  G 


3 


ax 


~7  d  F2  -3 

e  ^  "  B 


(  3.88  ) 
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By  equating  terms  of  like  coef f icients,  the  following 
steady-state  results  are  obtained 


a3  : 


r7 
G  : 


=  1 


(  3.89  ) 


^2 

ax  ax2ay 


ax 

ax2  ay 

axz 

1_  + 

dF  l  az 

F2 

i 

ay 

ax  ax 

2ay 

a3pi  »2 

ax3 

ay 

axz  ay 

azF2 

.  , 

2  ay 

ax2 

ax2 

(  3.90  ) 


=  0 


(  3.91  ) 


(  3.92  ) 


Eq  (  3.89  )  can  be  integrated  four  times  to  obtain  an 
expression  for  (X)  . 


X*  x3  x 2 

Fi(X)  =  +  ci  TT  +  C2  T~  +  c3  x  +  c4  • 


(  3.93  ) 


By  applying  the  boundary  conditions  (  Eqs  3.77  -  3.78  ),  the 
constants  of  integration  can  be  determined.  The  final  result  for 
VX)  is 


X  2 

F, (X)  =  2—  (  X  -  l  ) 

1  4 : 


(  3.94  ) 
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An  expression  -for  can  now  be  found  by  differentiating  Eq 
(  3.94  ),  twice,  with  respect  to  X  . 


d2F 


W. 


ar 


1  -  5  {  1(2  -  X  +  E  } 


(  3.95  ) 


From  Eq  (  3.90  ),  and  noting  that  aF^/aY  =  0  from  Eq 


(  3.94  ) ,  leads  to 


aV 


ax 


=  o 


(  3.96  ) 


Integrating  this  equation  four  times  gives  the  result  for 


F2<X)  . 


1  3 

F2<X>  "  61  X 


^  x2  +  C3X  +  C4  .  (  3.97  ) 


Applying  the  boundary  conditions,  the  constants  of 
integration  are  found  to  be: 


C,  =  C.  =  0 

3  4 


C,  =  -3  C, 


(  3.98  ) 


(  3.99  ) 


and 


C.  =  -2  C,  , 


(  3.100  ) 


which  can  only  be  true  if  Cj  =  C2  =  0  .  Therefore, 


F„(X)  =  0  , 


(  3.101  ) 


which  further  implies  that  W2<X)  =  0 
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With  these  results,  and  substituting  Eq  (  3.94  )  into 
(  3.86  )  and  Eq  (  3.95  )  into  (  3.87  ),  the  -final  form  of  the 
anaylitical  solution  for  the  energy,  stream  function,  :nd 
vorticity  are: 

T  =  1  -  X  (3. 102  ) 

2 

^3  X  7 

F  =  G  (  X  —  1  j  (3. 103  ) 

W=|(X2-X+^).  (  3.104  ) 

These  one  term  expansions  are  good  approximations  for 
modeling  pre-transi t i onal  flow  behavior  for  the  limiting 
condition  of  Pr  -*  0  .  In  their  present  form,  however,  one  could 
be  misled  to  conclude  that  they  are  applicable  throughout  the 
entire  vertical  slot.  This  is  not  true  because  the  original 
partial  differential  equations  (  Eqs  3.68  to  3.70  ),  from  which 
these  expressions  were  derived,  are  not  independent  of  Y  as  the 
analytical  solution  suggests.  Whereas,  the  analytical  solutions 
are.  The  best  comparision  with  the  numerical  solution  should 
occur  at  half  the  vertical  height  (  Y  =  1/2  ),  where  the  flow 
stream  and  vorticity  contours  for  a  single  cell  are,  to  a  good 
approximation,  constant  in  the  Y-direction. 

Boundary  Layer  Equations  for  Hi gh  Prandtl  Number  (  Pr  >  10.0  ) . 

Using  the  same  reasoning  as  for  the  low  Prandtl  number 
equations,  the  dimensionless  velocity  components,  vorticity. 
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stream  -function,  etc.  are 


u 

v 

W 

f 

T 

t 

G 


„  1/4  - 
Ra  u 


D  1/4  ~ 
Ra  v 


Ra3/4  W 


1/4  ~ 
Ra  -f 


-1/2  ~ 
Ra  t 


„  1/4  K 
Ra  G 


(  3.105  ) 

(  3. 106  ) 

(  3.107  ) 

(  3.108  ) 

(  3. 109  ) 

(  3.110  ) 

(  3.111  ) 


where,  again,  the  tilda  over  the  variable  signifies  a  scaled 
dimensionless  quantity.  By  substituting  Eqs  (  3.105  -  3.111  ) 


into  Eqs  (  2.44  -  2.46  ),  and  neglecting  the  terms  of  order 
~2  1/2  2  2 

G  /Ra  ,  a  /ay  ,  the  governing  equations  become 


Energy: 


B2  2- 

at 


r  /  51  -  *1  £L  \  -  9  1 

b  \  ax  av  ay  ax  j  a%2 


(  3.112  ) 


Vortici ty: 


G2 

at 


G 


f  w  aw 
\  ax  av 


»v»  ~  ^ 

aw  £f  1 
ax  ay  ) 


=  Pr 


+ 


Sh  \ 
»x2  J 


(  3.113  > 
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Stream  Function 


1 


W 


(  3.114  ) 


Note  that  the  vorticity  equation  (  Eq  3.113  )  is  now 
diffusion-dominated  -for  high  Prandtl  numbers.  The  boundary  and 
initial  conditions  remain  the  same. 

Boundary: 


f (O,  Y)  =  f ( 1 , Y)  =  f<X,0)  =  f(X,l>  =  0 


ax 


x=o,  1 


av 


=  o 


Y=0 ,  1 


T  (0,Y)  =  1,  T ( 1 , Y)  =  0 


(  3.115  ) 

(  3.116  ) 

(  3.117  ) 


aj 

ay 


=  o 

Y=0, 1 


(  3.118  > 


Initial : 


X=0, 1 


G2 


W  . 


i  =  W  =  T  =  0  , 


except  at  the  left  vertical  wall  where 


T 


1 


(  3.119  ) 


(  3. 120  ) 


(  3.121  ) 


Limiting  Equations  for  Infinite  Prandtl  Number .  To  make 
the  study  complete,  one  would  also  like  to  investigate  the 


III 


20 


behavior  o-f  -fluids  as  the  Prandtl  number  becomes  in-finite.  To 
accomplish  this,  the  high  Prandtl  number  equations  are  simpli¬ 
fied  for  the  limiting  condition  of  Pr  -»  »  . 


First,  let  Pr  -*  oo  and  assume  that  temperature  remains  0(1). 
Now,  assume  the  vorticity  and  stream  functions  have  the  form 


W  =  7>i  Wi  + 


f  =  r)2  ^ 


T  =  T  + 


(  3.122  ) 

(  3.123  ) 

(  3.124  ) 


Substituting  these  expressions  into  Eqs  (  2.44  -  2.46  ),  the 
energy,  vorticity,  and  stream  function  equations  become 


Energy: 


2  - 


a  «T  ^  a  f  af  91  «T  1  _ 

G  ~  ^2  6  \  ax  ay  aY  ax  j 


a2T 

ax2 


(  3.125  ) 


Vorticity: 


x2  aw 

6  — 

A  at 


p  /  af  aw  _  aw  af  1 
77 1  v2  b  \  ax  ay  ax  ay  J 


=  Pr 


f  ?  ar  ,  a2w  \ 

1 G »  +  j? ) 


(  3.126  ) 


Stream  Function: 


a2f 


v?  — ~  =  r).  e2  w 
err  1 


But,  from  the  stream  function  equation,  Eq  (  3.114  ), 


(  3. 127  ) 
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(  3.128  ) 


Therefore,  in  the  limit  as  pr  ■*  a>  ,  ~  ~  0(1)  .  Also, 

for  this  limit,  the  vorticity  equation  reduces  to  the  Stokes- 
like  momentum  equation.  The  energy,  vorticity,  and  stream 
function  equations  simply  become 


Energy: 


G2  *1 

at 

Vorticity: 


/>  fv  IV  V  s  7 IV 

~fdf£T  _  «  3T  1  £T 

|  «?X  <?Y  dY  dX  J  ax2 


(  3.129  ) 


G 


dT 

dx 


o 


(  3.130  ) 


Stream  Function: 


(  3.131  ) 


The  boundary  and  initial  conditions  are  similar  to  those  in  Eqs 
(  3.115  -3. 121  ) . 


Steady-State  Analytical  Solution  for  Pr  -»  oo  .  Using  the 
same  procedure  as  before,  the  steady-state  perturbation  solution 
for  vorticity,  stream  function,  and  temperature  are  found. 
Vorticity: 

W  =  G  W1  (3. 132  ) 

Stream  Function: 

f  =  G  3  f  (3. 133  ) 
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Temperature 


where 


A 


1 


T  =  (  1  -  X  )  +  G8  ( 

5>  -  HxZ-x*H 

-  1  2  2 
*1  =  I4  X  (  X  -  1  > 

-1  A2 
T  ss  _ _  A  +  _ 

2  96  1  34560 

-X11  X10  -X9 

415800  +  75600  +  45360 

-x5.  X4 
+  6900*  15120  + 


+  10800 

-X 

387701 


(  3-134  ) 

(  3.135  ) 

(  3.136  ) 

(  3.137  ) 

(  3.138  ) 


_  —4  y 1 1  2  Y10  .  -5  y9  3  y8  -5  Y7 

2  55  5*  6*  4  *  2T  * 


1 2  x5 

10  * 


(  3.139  ) 


The  boundary-layer  equations  developed  in  this  chapter  -for 
low  Prandtl  number,  Pr  ■*  0,  high  Prandtl  number,  and  Pr  +  ao 
were  solved  numerically  to  obtain  stream  -function,  vorticity, 
and  temperature  contours  -for  the  vertical  slot  geometry.  The 
numerical  techniques  employed  are  discussed  in  the  Chapter  IV. 
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IV.  Numerical  Anal ysi s 

There  are  two  topics  of  discussion  in  this  chapter.  The 
■first  describes  the  -f inite-di-f -f erencing  method  used  in  solving 
the  coupled  set  of  governing  partial  differential  equations. 
These  equations  are  the  2-D,  unsteady,  boundary-1 ayer  equations 
developed  in  Chapter  III.  The  second  topic  deals  with  the 
computational  procedure  employed  to  solve  the  resulting 
f ini te-dif f erence  equations.  Here,  the  iteration  sequence, 
convergence  criterion,  and  relaxation  parameters  for  the 
numerical  procedure  are  described. 

The  basic  numerical  approach  was  adopted  from  Fant 
(  1987  ).  In  this  study,  however,  the  geometry  was  changed  from 
that  of  an  annulus  to  a  vertical  slot.  The  governing  equations 
were  converted  from  polar  to  Cartesian  form. 

In  computing  the  streamwise  (  Y-direction  )  nonlinear 
convective  terms,  a  corrected,  second-order  upwind  scheme  was 
used  for  the  boundary-1 ayer  equations.  For  the  X-coordinate 
convective  terms,  a  corrected,  second-order  central-dif f erencing 
scheme  was  used.  This  approach  ensures  numerical  stability  when 
solving  highly  convective  flow  problems.  The  unsteady  form  of 
the  equations  provided  the  opportunity  to  capture  both  steady 
and  unsteady  flow  behavior  when  solved  in  a  time  accurate 
fashion. 
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Numerical  Method 


The  2-D,  unsteady,  boundary-1 ayer  equations  presented  in 
Eqs  (  2.30  -  2.32  )  and  (  2.44  -  2.46  )  are  discretized  using 
f ini te-dif f erencing  techniques.  In  this  -form,  the  governing 
equations  are  well  suited  for  solving  on  a  high-speed  computer 
system.  The  computational  domain  consists  of  a  rectangular  mesh 
with  more  nodes  in  the  vertical  (  Y  )  than  in  the  horizontal 
(  X  )  direction  (  Figure  2  ).  Non-di mensi onal i z i ng  the  X  and  Y 
spacial  coordinates  with  L  and  1,  respectively,  transforms  the 
geometrical  domain  into  a  square  box  with  all  sides  equal  to 
1.0.  The  computer  code  used  to  perform  the  computations  had  the 
capability  to  use  variable  increments  in  the  computational 
domain.  Although  the  geometry  is  symmetrical,  vertical  and 
horizontal  symmetry  were  not  assumed  in  the  calculations. 

Finite-Difference  Expressions  for  Variable  Increments.  In 
the  governing  equations,  Taylor  series  expansions  are  used  to 
transform  the  derivatives  into  finite-difference  expressions. 
Standard  forward-,  backward-,  and  central -difference  expressions, 
including  formal  truncation  c^^ors  for  variable  increments,  are 
given  in  Fant  (1987). 
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The  unsteady  derivatives  are  approximated  with  a  stable. 


-f orward-dif -f erence  molecule 


d<p 

at 


o 


At 


0 (At )  , 


(4.1) 


where  n  is  the  time  level  at  which  the  variable  <p  is  evaluated. 


Finite-Difference  Equations  -for  the  Dependent  Variables. 

All  spacial  derivatives  in  the  governing  equations  are  second- 
order,  centrally-  (  or  upwind-  )  differenced.  This  includes  the 
convective  terms  in  the  vorticity  and  energy  equations  which  are 
represented  by  a  first-order,  upwind  expression  together  with  a 
correction  term  for  second-order  accuracy.  The  second-order 
terms  are  split  in  this  manner  in  order  to  enhance  numerical 
stability  when  resolving  secondary  flow  behavior.  This 
procedure  is  explained  more  fully  in  Fant  (  1987  ).  For  the 
limiting  case  where  Pr and tl  number  goes  to  zero,  the  energy 
equation  (  Eq  3.68  )  decouples  from  the  vorticity  and  stream 
function  equations  <  Eqs  3.69  and  3.70  ).  Therefore,  the  energy 
equation  can  be  solved  directly  while  the  other  two  have  to  be 
solved  numerically  in  a  coupled  manner. 

2-D  Navi er-Stokes  Equations.  The  nonlinear  convective 
terms  in  the  Navi er-Stokes  equations  can  be  expressed  as 


) 
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2  x  ax 


=  (  X 


x=o 


-  N  .  P 


<^1  -  H  03  "  <  1  +  H  }  <^0 
hb  +  h* 


} 


+  (  X 


IM  >  p 


4>,  ~  4>f 
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(  4.2  ) 
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*o  ~  *4 
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kb  +  kf 


} 
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-} 


(  4.3  ) 


where  tp  represents  the  dependent  variables  W,  -f,or  T,  and 


hf  kf 

H  =  —  and  k  =  —  . 


<  4.4  ) 
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Eqs  (  4.2  )  and  (  4.3  )  are  a  combination  of  first-order, 
upwind-difference  components  with  the  added  corrrections  to 
bring  the  differencing  up  to  second-order  accuracy.  Using  this 
differencing  technique,  the  energy,  vorticity,  and  stream 
function  boundary-layer  equations  developed  in  Chapter  III  can 
be  manipulated  into  the  following  general  form: 


A2 


{ 


~  2 
A  G 


dtp 

dt 


d2*  \ 
dX2  ) 


x  ** 
x  dx 


2  u  —t- 


d<P 

d'i 


s  = 


(  4.5  ) 


The  coefficients  +  |\|  and  +  |^|  were  added  to  ensure  that 
stable  differencing  ' into-the-wind '  is  enforced.  G,  t,  X, 
A2,  tpn  and  A  have  different  meanings  depending  upon  the 
application  of  the  boundary-1 ayer  equation  (  high,  low,  or 
limiting  Pr  ) . 


Low  Prandtl  Number  <  Pr  <0.1  ) : 


G  =  G,  t  =  t 


2  \  = 


2  — 


Ai  *£ 

A1  dv 


-A1  *1 

dx 


(  4.6  ) 


(  4.7  ) 


(  4.8  ) 
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A1  = 


Pr  G,  <p  =  T 

Pr  G,  <p  *  W 

0,  <f>  *  ? 


A2  = 


f1’ 

4  Pr* 


0  =  T 
<*>  =  W 


1/  G  ,  0  =  i 


(  4.9  ) 


A  = 
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^  0,  <*>  =  ? 
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r  o* 

p  dT 
_G  3X’ 
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*  =  T 
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(  4.10  > 


Zero  Prandtl  Number  (  Pr  -*  0  )  : 


G  =  6,  t  =  t 


(  4.11  ) 
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(  4.12  ) 


(  4.13  ) 


(  4.14  ) 


O, 

1. 

0, 


<P  «  T 
<P  =  W 
*  =  F 


S  = 


O,  4>  -  T 

-  G,  <fi  =  W 

-  W,  <p  =  F 


(  4.15  ) 


High  Prandtl  Number  (  Pr  >  10.0  ): 

G  =  6,  t  =  t 


(  4.16  ) 
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With  the  above  results,  the  boundary-1 ayer  equations  can 


be  written  in  the  following  finite-difference  form: 

2 


C  4> 

o  o 


n  +  1 
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(  4.26  ) 


with  the  variable  <p  denoting  the  boundary— 1  ayer  dependent 
variables  T,  W,  F,  and 


A2  A  G  2  (A2) 
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(  4.28  ) 


(  4.29  ) 


C,  =  2  (A2) 


{  V  h!  *  V  } 
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(  4.30  ) 
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(  4.31  ) 
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where 


K, 


'f  1 


Kl 


and 


A  =  k. 


bl 


bl 


B  =  k, 


'f  1 


Boundary  Conditions.  For  the  vertical  slot,  the 
difference  form  of  the  boundary  conditions  are 

Energy: 


Tn+1  ,  .  Tn+1  - 

T .  =  1  and  Tk._  =  0 

l , j  NR , j 


31 

ay 


i  ,  1 
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ay 
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i  ,NS 


Vortici ty: 
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n+1 
1 ,  J 


2  F 


n+1 

2,J 


(  G  hj  ) 
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NR,  j 


2  F 
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1,1  i , NS 


(  4.34  ) 

(  4.35  ) 

f i ni te— 

(  4.36  ) 

(  4.37  ) 

(  4.38  ) 

(  4.39  ) 

(  4.40  ) 
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Stream  Functions 


n  +  1  _  p.n+1  _  _n+l  _  _n+l 

l,j  NR,j  i,l  i,NS 


(  4.41  ) 


where  i  goes  from  1  to  NR,  and  j  from  1  to  NS  within  the 
geometric  boundaries.  Additional  nodes  exist  outside  these 
boundaries  (i=0,j=0,i=  NR  +  1,  and  j  =  NS  +  1)  for 
computational  purposes. 


It  should  be  noted  that  the  computer  code  used  to  solve  the 
vertical  slot  problem  was  originally  developed  for  an  annulus. 
Since  an  annulus  has  only  two  boundaries,  inner  and  outer,  a  set 
of  finite-difference  equations  was  needed  for  the  adiabatic 
horizontal  walls.  As  mentioned  previously,  the  finite- 
difference  expressions  for  the  dependent  variables  are  second- 
order  accurate.  This  suggests  the  use  of  a  second-order 
polynomial  to  represent  the  adiabatic  boundary  conditions  along 
the  top  and  bottom  horizontal  walls.  Using  a  Taylor  series 
expansion,  a  second-order  accurate,  finite-difference  expression 
was  obtained. 


Along  the  bottom  horizontal  (  j  =  1  >  wall. 
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i  ,  3 
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=  0  (  4.42  ) 
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from  which  the  temperature  is  simply 


4  T 


i,2 


i  » 3 


i  ,1 


3 


<  4.43  ) 


For  the  top  horizontal  wall  (  j  =  NS  ) 


dl 

ay 


i  ,  NS 


Ti,NS-2  4  Ti ,NS-1  +  3  Ti ,NS 
2  AY 


=  0 


(  4.44  ) 


and , 


4  T 

i ,NS-1 


-  T 


i , NS— 2 


i  ,  NS  3 

To  satisify  the  adiabatic  constraints,  the  fallowing 
conditions  also  apply: 


Ti , 2  Ti,0  and  Ti ,  NS+1  Ti ,NS— 1 


(  4.45  ) 


(  4.46  ) 


The  boundary  conditions  along  the  left  vertical  wall  were 
obtained  by  using  a  Taylor  series  expansion  and  then  solving  for 
the  second  derivative. 


=  F. 
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(  4.47  ) 


From  the  no-slip  constraint,  F^  and  aF^/  ax 


=  0  .  Therefore, 


with  AX  =  h.  , 


aV 


ax‘ 
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<  ) 


(  4.48  ) 
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and 


W  = 


2  F, 


(  G  hj  > 


(  4.49  ) 


The  right  vertical  wall  boundary  condition  was  handled  in  a 
similar  manner. 


The  vorticity  boundary  conditions  for  the  horizontal  walls 


are 


ln+l 

i  i 


W 


n+1 
i  ,  NS 


0  . 


(  4.50  ) 


These  conditions  are  a  result  of  the  boundary-layer  approxi¬ 
mation,  where 


d2F  J  f  £  ] 

*X2  ~  t  ax  J 


(  4.51  ) 


But,  at  the  walls,  the  no-slip  constraint  requires  that 

-3^-  =  0.  This  implies  that  W  must  also  be  zero.  However,  to 

ox 

match  the  physics  of  the  problem,  as  well  as  ensure  that  W  =  0 
at  the  top  and  bottom  walls,  requires  that 


ln+1 
i  » 


-  W 


n+1 

i  ,2 


and 
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n+1 

i , NS+1 


-  W 


n+1 
NS— 1 


(  4.52  ) 


Computational  Procedure. 

The  system  of  coupled,  finite-difference  equations 
described  in  the  preceding  section  were  solved  implicitly  in 
time  using  a  point-i terati ve,  Gauss-Siedel  scheme  with 
under-rel axation.  A  computer  code  was  developed  for  each 
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Prandtl  number  range  investigated:  low  Prandtl  number,  Pr  -►  0, 
Prandtl  number  of  order  one,  and  high  Prandtl  number.  In  each 
case,  the  dependent  variables  were  solved  numerically  by 
repeated  iterations  of  the  governing  equations  at  a  given  time 
1 evel . 


At  the  beginning  of  each  investigation,  the  initial 
conditions  for  temperature,  vorticity,  and  stream  function  were 
set  to  zero  everywhere  within  the  slot,  except  for  the  iso¬ 
thermal  condition  along  the  left  vertical  wall.  In  addition, 

~  '*• 

the  scaled  gap  parameter  (  G  or  G  )  was  set  to  1.0  and  a  node 
generating  program  was  employed  to  set-up  the  computational 
domain.  The  output  of  the  initial  computer  run;  consisting  of 
temperature,  vorticity,  and  stream  function  distributions,  was 
then  used  as  initial  conditions  for  the  next  run.  In  this  new 
run,  the  scaled  gap  was  increased  to  a  slightly  higher  value  and 
the  procedure  repeated  until  transition  to  multicells  occured. 

Iteration  Sequence  and  Convergence  Cri teria.  The  finite- 
difference  form  of  the  boundary-1 ayer  equations  were  iterated 
in  the  following  order: 

i .  energy 

ii.  vorticity 

iii.  stream  function 

except  for  the  small  Prandtl  number  limiting  equations  (  Pr  -»  0  )  , 
where  an  analytical  expression  was  obtained  for  energy. 
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This  sequence  was  repeated  until  successive  iterations  were 
within  a  prescribed  tolerance.  The  relative  convergence 
constraint  -for  the  maximum  modulus  of  the  difference  between 
energy,  vorticity,  and  stream  function  is  defined  as 

<  1  x  10_A  (  4.53  ) 

where,  again,  <p  represents  temperature,  vorticity,  or  stream 
function.  When  Eq  (  4.53  )  was  satisified,  the  numerical 
solution  was  considered  converged.  In  this  equation,  m  denotes 
the  iteration  level  and  n  denotes  the  time  level. 


Max 


,m+l,n+l  .m,n+l 

<P  <P 


<  m+1 ,n+l  } 


Relaxation  Parameters.  To  control  the  rate  of  convergence, 
two  independent  relaxation  parameters,  and  ,  were 

incorporated  into  the  finite-difference  expressions.  The 
parameter,  ,  was  employed  in  conjunction  with  the  vorticity 

boundary  conditions,  namely 

„p.m,n+l 
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along  the  right  wall  where  h  represents  the  linear  spacing 
between  nodes.  The  values  used  to  stabilize  the  numerical 
computations  were  generally  in  the  range  . 15  <  0^  <  .3  .  The 

second  parameter,  was  used  with  the  second-order  upwind 

correction  terms  which  were  based  on  the  new  time  level,  n+1  . 
The  numerical  values  for  vorticity,  at  each  interior  node,  were 
computed  according  to 


w<n+l  ,n+l 
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,  Em,n+1 


} 


(  4.56  ) 


The  correction  term,  E,  becomes  a  weighted  average  of 
Gauss-Si edel  (  at  the  current  level  )  and  Jacobi  (  at  the 
previous  level  )  iterations  for  values  of  <  1.0  .  Values  of 
. 15  <  02  <  .3  were  generally  employed  in  the  computations.  For 
further  details  on  this  procedure,  see  Fant  (  1987  ). 

Numerical  Detai 1 s  and  Pi scussi on.  The  majority  of  the 
computational  runs  were  made  using  31  horizontal  and  62  vertical 
nodes.  Because  of  the  large  gradients  expected  near  the  left 
and  right  vertical  boundaries,  the  number  of  nodes  were  slightly 
concentrated  in  those  regions.  This  helped  to  improve  numerical 
accuracy. 
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A  -few  cases  were  run  with  a  31  x  102  mesh.  This  was  done 
■for  two  reasons:  -first,  to  determine  the  influence  of  mesh  size 
on  triggering  the  onset  of  cell  formation;  and  second,  to 
concentrate  more  nodes  in  the  center  of  the  box  to  resolve  the 
temperature  inversions  expected  for  the  high  Prandtl  number 
i nvest i gat i ons. 

As  explained  earlier,  each  investigation  was  started  by 
computing  the  flow  field  for  G  =  1.0  .  The  output  from  each 
computational  run  served  as  the  input  for  the  next  with  G 
incremented  to  a  new  value.  The  procedure  was  repeated  until 
the  formation  of  multiple  cells  appeared  in  the  results. 

A  scaled  time  step.  At,  was  used  in  each  Prandtl  number 

~  1/2 

investigation.  Since  At  =  <  Ra/Pr  )  At  ,  the  actual  time 
step  is  much  smaller.  Generally,  a  At  of  0.5  was  used  and  the 
total  number  of  time  steps,  r,  was  about  500  . 

From  earlier  studies  done  by  Fant  <  1987  )  and  Bennett 
(  1988  ),  a  suitable  range  of  values  for  G,  where  cell  formation 
might  occur,  was  surmized.  Initially,  G  was  incremented  by  as 
large  a  value  as  possible  for  each  succeeding  run.  This 
increment  generally  ranged  from  1.0  <  G  <  3.0  .  Once  cell 
formation  developed,  G  was  reset  to  the  previous  value  where  a 
single  cell  was  still  present.  The  results  from  the  previous 
(  single  cell  )  computational  run  were  then  used  as  input,  but 
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with  the  increment  for  G  reduced  <  AG  ~  . 1,  .2,  etc.  ). 
Additional  runs  were  then  accomplished  until  mul ticel lular 
formation  ocurred  again,  with  the  value  for  §  now  more 
accurately  determined. 

Prior  to  the  onset  of  multiple  cell  formation,  the  CPU  time 
required  to  obtain  the  results  for  each  run  was  usually  on  the 
order  of  30  minutes  to  an  hour.  However,  when  the  flow  began  to 
transition,  the  number  of  iterations  for  convergence  and  the  CPU 
time  increased  significantly.  This  was  particularly  true  for 
the  high  Prandtl  number  cases  where  CPU  time  was  approximately 
24  hours  per  run  for  Pr  =  1000  . 

The  validity  of  this  last  statement  is  suspect,  however, 
since  major  modifications  were  made  to  the  VAX  computer  system 
which  may  have  slowed  the  processing  speed  considerably.  These 
modifications  were  done  while  the  study  was  in  progress.  Prior 
to  this,  the  computational  runs  proceeded  at  a  fairly  rapid  pace 
(  several  accomplished  each  day  ).  But,  after  that,  things 
slowed  down  dramatical ly.  One  run  required  over  90  hours 
elapsed  (  24  hours  CPU  )  time  to  complete.  Thus,  the  amount  of 
studies  accomplished  for  the  high  Prandtl  number  regime  were 
reduced. 
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V 


Results  and  Discussion 


In  this  chapter,  the  various  flow  fields  computed  for  the 
low,  high,  and  limiting  Prandtl  number  regimes  are  examined. 

This  includes  temperature,  vorticity,  and  stream  function 
distributions  throughout  the  vertical  slot.  In  each  case,  the 
hot,  isothermal  wall  appears  on  the  left,  while  the  right 
vertical  wall  is  maintained  at  zero.  To  establish  validity  in 
the  numerical  procedure,  the  results  are  compared  to  analytical 
solutions  obtained  for  steady-state  conditions. 

Zero  Prandtl  Number. 

Needless  to  say,  a  fluid  with  a  Prandtl  number  equal  to 
zero  does  not  exist  in  the  real  physical  world.  It  is  a 
hypothetical  fluid  whose  thermal  diffusivity  approaches 
infinity.  In  this  case,  the  energy  equation  reduces  to  the 
simple  steady-state  conduction  equation.  For  the  vertical  slot, 
this  is  just  a  linear  temperature  distribution  across  the  width 
of  the  slot  for  the  boundary  conditions  imposed. 

Since  the  thermal  diffusivity  is  infinite,  thermal 
perturbati ons  cannot  exist  because  they  are  instantaneously 
diffused  throughout  the  slot.  Thus,  any  instabilities  generated 
are  purely  hydrodynamic  in  nature.  The  way  in  which  the  flow 

A 

develops  for  increasing  G  is  shown  in  the  stream  function  plots. 
Figure  3  (  a  -  e  ) . 
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Figure  3a  shows  the  streamlines  for  the  steady-state 
condition  where  G  =  3.0  .  The  motion  of  the  fluid  is  always  in 
a  clockwise  manner  throughout  the  slot.  Notice  that  even  when 
the  flow  is  stable,  weak  convective  cells  have  formed  at  the 
extreme  top  and  bottom  of  the  cavity.  Cells  form  in  this  region 
first  because  the  boundaries  force  the  fluid  to  turn.  Since 
momentum  diffusivity  is  so  low,  the  motion  does  not  dissipate 
quickly  enough  causing  some  of  the  fluid  particles  to  continue 
in  a  circular  path.  Lee  and  Korpela  (  1983  )  obtained  a  similar 
one  cell  result  using  the  complete  Navier — Stokes  equations. 
However,  as  mentioned  by  the  authors,  their  plot  did  not  show 
the  weak  cells  at  the  top  and  bottom  because  of  insufficient 
resolution.  As  G  increases,  these  cells  become  stronger 
(  Figure  3b  and  c  )  where  forced  turning  of  the  fluid  by  the 
physical  boundary  continues  to  contribute  to  their  development. 

Finally,  at  G  -  4.0  (  Figure  3d  >,  transition  begins  to 
take  place  and  multiple  cells  begin  to  form.  Cell  formation 
takes  place  near  the  top  and  bottom,  first,  and  then  toward  the 
center.  Only  a  moderate  increase  in  G,  from  4.0  to  4.2,  was 
needed  to  trigger  complete  multicellular  formation  (  Figure 
3e  ) .  A  total  of  nine  cells  are  present  with  the  strongest 
cells  appearing  near  the  horizontal  walls  and  the  weaker  ones  in 
the  center.  Figure  3  also  shows  that  the  entire  fluid  flow  is 
slightly  skewed  toward  the  upper  right  cold  wall  and  toward  the 
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lower  left  hot  wall,  respectively. 


Figures  4  and  5  show  the  comparison  of  the  numerical 
results  with  the  analytical  solution  (  Eq  3.103  )  obtained  for 
the  stream  function,  at  the  horizontal  centerline,  for  the 
steady-state  condition  of  G  =  3.0  .  The  agreement  between 
numerical  and  analytical  results  is  excellent.  Figure  5  is  more 
revealing  in  that  it  shows  that  the  largest  relative  error 
(  ~  5  7.  )  between  the  numerical  and  analytical  solutions  occurs 
very  close  to  the  ve  fical  boundaries.  This  is  acceptable, 

because  it  represents  the  difference  between  two  numbers  very 

-4 

close  to  zero  (  ~  10  ). 

Equally  as  interesting,  is  the  development  of  vorticity 
throughout  the  cavity.  This  is  shown  in  Figures  6  (  a  -  e  ) , 
which  are  contour  plots  of  constant  vorticity.  Vorticity  is  in 
the  counter — clockwise  direction  (  positive  )  along  the  vertical 
walls,  and  clockwise  (  negative  )  in  the  cental  region  of  the 
slot 


Note  how  the  vorticity  contours  arrange  themselves  into  a 
pattern  similar  to  that  for  the  stream  function.  Notice,  also, 
that  the  vorticity  tends  to  concentrate  in  specific  nodal 
regions  along  the  vertical  walls  (  Figure  6e  ).  These 
concentrations  form  near  the  ends  of  the  walls,  first,  and  then 
towards  the  center  (  Figure  6c,  d,  and  e  ).  Also,  the 
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concentrations  along  the  left  wall  do  not  align  with  their 
counterparts  on  the  right.  This  may  be  a  result  of  the  skewed 
•fluid  motion  mentioned  earlier.  Lastly,  the  contours  that 
intersect  the  horizontal  walls  represent  zero  vorticity.  This 
is  the  point  where  vorticity  changes  direction,  from  counter — 
clockwise  to  clockwise  or  vice-versa. 


It  is  interesting  to  note  that  there  are  no  vorticity 
concentrati ons  along  the  horizontal  walls.  This  may  be  due  to 
the  boundary-layer  approximation  to  the  Navier — Stokes  equations 
as  Ra  -*  oo  and  G  -»  0  where 


G2  a2w 
Ra1/2  *Y2 

and  the  associated  boundary  conditions  are  neglected. 


Comparisons  between  analytical  (  Eq  3.104  )  and  numerical 
results  at  the  horizontal  centerline  are  shown  in  Figures  7  and 
8  .  Again,  the  comparison  is  quite  good.  However,  this  time 
the  maximum  error  occurs  not  near  the  vertical  wall,  but  at  the 
region  where  the  vorticity  changes  direction  (  Figure  8  ). 
Again,  the  error  is  small  (  ~  3.0  V.  )  and  represents  the  error 
between  two  numbers  very  close  to  zero. 


As  expected  for  the  zero  Prandtl  number  case,  the  isotherm 
contours  were  simply  straight  vertical  lines  and  are  of  no 
interest.  Therefore,  they  are  omitted  from  this  study. 
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Low  Prandtl  Number ■ 

Li  quid  Metal s  (  Pr  ~  0.02  ) .  When  the  Prandtl  number  of 
the  -fluid  is  increased  to  that  o-f  liquid  mercury,  very  little 
changes.  Figure  9  shows  the  stream  -function  contours  -for  a 
fluid  o-f  Pr  =  0.02  .  The  same  pattern  of  cell  development  as 
that  for  Pr  =  0  was  observed  with  cell  formation  occuring  near 
the  horizontal  walls,  first,  and  then  in  the  central  region. 

Figure  9a  shows  the  fluid  flow  for  G  =  1.0  .  Even  for  this 
low  value,  the  single  cell  shows  a  definite  necking  down  region 
near  the  horizontal  walls  where  weak  convective  cells  are 
begining  to  form.  For  G  =  1.5,  multicellular  development  is 
well  in  progress  (  Figure  9b  )  and  is  complete  when  G  =  1.6 
(  Figure  9c  ) . 

nj  A  1/4 

Recalling  that  G  =  G  Pr  ,  and  substituting  the  value  for 
G  at  which  multicellular  formation  occured  for  the  zero  Prandtl 
number  case  (  ~  4.2  ),  leads  to  G  =  1.57  which  is  in  good 
agreement  with  the  actual  value  obtained.  Further  increases  in 
G  did  not  cause  an  increase  in  the  number  cells  that  formed,  but 
only  increased  their  strength.  As  a  result,  steady-state 
conditions  were  maintained  over  the  range  of  G  that  was  studied. 
Notice,  also,  that  nine  cells  formed  with  the  strongest  cells 
near  the  horizontal  boundaries;  just  as  in  the  case  for  Pr  =  0  . 
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Figures  10  and  11  show  the  comparison  between  numerical  and 
analytical  (  Eq  3.40  )  results.  Again,  the  comparisons  are  good 
and  show  similar  trends. 


Vorticity  also  developed  in  the  same  manner.  This  is  shown 
in  Figure  12.  Figure  12e  clearly  illustrates  the  tendency  for 
vorticity  to  concentrate  in  specific  regions  along  the  vertical 
walls  and  has  trends  similar  to  the  Pr  =  0  case.  And  as  before, 
the  analytical  (  Eq  3.41  )  and  numerical  comparisons  for 
vorticity  were  very  g  _<od  (  Figures  13  and  14  ). 


Temperature  contours,  however,  were  no  longer  vertical 
straight  lines  <  Figure  15  ).  For  G  =  1.8  (  Figure  15e  ),  the 
contours  showed  a  distinct  oscillation  in  regions  where  cell 
formation  took  place.  This  is  because,  with  a  finite  Prandtl 
number,  thermal  diffusion  is  no  longer  instantaneous.  There¬ 
fore,  the  temperature  distribution  is  slightly  influenced  by  the 
convective  motion  of  the  fluid.  Similar  results  were  reported 
by  Lee  and  Korpela  (  1983  >. 


As  expected,  the  analytical  (  Eq  3.42  )  and  numerical 
solutions  for  the  single  cell,  steady-state  condition  compared 
very  well  (  Figures  16  and  17  ). 
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Prandtl  Number  o-f  Order  One. 


Ai r  (  Pr  ~  0.70  ).  With  air,  an  interesting  phenomenon  was 
observed.  For  G  =  3.9,  two  tiny  cells  -formed  in  the  center  o-f 
the  vertical  slot  (  Figure  18b  ).  Increasing  6  to  4.0  (  Figure 
18c  )  caused  an  immediate  formation  of  multiple  cells,  with  the 
stronger  cells  occurring  in  the  center  of  the  slot.  This  is  in 
direct  contrast  to  the  very  low  Prandtl  numbers,  where  the  cells 
formed  near  the  horizontal  boundaries,  first. 

For  air,  momentum  diffusivity  is  higher.  Therefore,  fluid 
motion  dissipates  more  quickly  and  forced  turning  of  the  fluid 
by  the  boundaries  no  longer  has  the  same  effect.  As  a  result, 
cells  no  longer  form  adjacent  to  the  horizontal  boundaries, 
reducing  the  total  number  of  cells  from  nine  to  seven.  Further¬ 
more,  cells  form  in  the  center  of  the  slot,  first,  due  to  shear¬ 
ing  caused  by  the  opposing  movment  of  adjacent  fluid  layers. 

A  total  of  seven  distinct  cells  are  visible  in  Figure  18c. 
However,  as  G  was  further  increased,  the  lower  cell  disappeared 
while  the  remaining  ones  became  stronger  (  Figures  18d  and  e  ). 
Again,  these  results  are  very  similar  to  that  obtained  by  Lee 
and  Korpela  (  1983  )  using  the  complete  Navier — Stokes  equations. 
In  L'-'eir  study,  they  reported  that  as  the  cells  became  stronger, 
the  distance  between  them  increased.  The  same  happened  in  this 
study,  causing  the  cell  at  the  bottom  to  be  absorbed.  Notice 
that  this  time  however,  the  strongest  cells  are  not  in  the 
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regions  near  the  horizontal  boundaries,  as  was  the  case  for 
Pr  =  0  and  Pr  =  0.02,  but  are  located  in  the  center  of  the 
vertical  slot  where  shearing  of  the  fluid  is  strongest. 

Figure  19  shows  the  vorticity  distribution  for  air.  It  is 
interesting  to  note  that  the  contours  of  constant  vorticity 
verify  that  the  distance  between  cells  is  increasing.  Figure 
19c  clearly  shows  a  line  of  constant  vorticity  surrounding  the 
seven  convection  cells,  shown  in  Figure  18c  (  stream  function  ). 
Figure  19d  indicates  that  regions  of  constant  vorticity  have 
formed  between  the  top,  center,  and  bottom  cell  groups.  The 
only  way  for  this  to  happen  is  for  the  cells  to  be  spacially 
separating. 

Constant  temperature  (  Figure  20  )  contours  obtained  for 
air  were  very  similar  to  those  reported  by  Lee  and  Korpela 
(  1983  ).  Temperature  contours  tend  to  bend  toward  the  cold 
wall  along  the  upper  horizontal  surface,  away  from  it  along 
the  lower  wall,  and  also  oscillate  in  the  region  where 
convective  cells  form  (  Figure  20c,  d,  and  e  ).  The  only 
notable  difference  was  in  the  slight  temperature  reversal,  near 
the  horizontal  boundaries,  that  they  obtained. 

As  has  been  the  case  so  far,  analytical  (  Eqs  3.40  -  3.42  ) 
and  numerical  comparisons  for  the  stream  function,  vorticity, 
and  temperature  at  the  centerline  are  quite  good  and  also  show 
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the  same  trends  as  before  (  Figures  22-26  ). 

In  Chapter  III,  it  was  mentioned  that  because  the 
analytical  equations  are  a  -function  of  X  only,  one  could  be 
misled  to  conclude  that  they  are  valid  throughout  the  vertical 
slot.  To  verify  this  statement,  the  analytical  solution  was 
compared  to  the  numerical  solution  at  other  locations  a  distance 
of  20  percent  from  the  top  and  bottom,  respectively.  This 
corresponds  to  nodes  i,  12  and  i,  50  .  Figures  (  36  -  53  ),  in 
appendix  A,  show  that  the  analytical  solutions  for  temperature, 
vorticity,  and  stream  function  are  still  a  very  good  approxi¬ 
mation  for  distances  up  to  20  percent  from  the  horizontal  walls. 

Generally,  the  difference  between  the  numerical  and 
analytical  solutions  for  stream  function  and  temperature  was 
less  than  ten  percent.  Vorticity,  however,  was  the  exception 
with  an  error  of  20  percent.  Notice  that  in  the  lower  portion 
of  the  vertical  slot,  the  maximum  error  occurs  near  the  left 
boundary;  and  then  shifts  to  the  right  side  near  the  top  of  the 
slot  (  Figures  36-53  ).  The  reason  for  this  can  be  seen  from 
the  stream  function,  vorticity,  and  temperature  contours 
(  Figures  18-20  ).  The  analytical  solution  is  only  valid  at 
the  horizontal  centerline  of  the  slot  while  the  contours  for 
stream  function,  vorticity,  and  temperature  bend  toward  the 
lower  left  and  upper  right  vertical  walls,  respectively. 
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High  Prandtl  Number . 

Oils  (  Pr  ~  1000  )  .  For  high  Prandtl  numbers,  the  -flow  is 
momentum-diffusion  dominated  and  instabilities  that  result  are 
more  thermal  in  nature.  Bas.  -flow  heat  transfer  is  convection 
dominated  to  such  an  extent  that  distinct  boundary  layers  form 
on  each  wall  in  a  vertical  slot.  de  Vahl  Davis  and  Mallinson 
(  1975  )  and  Lee  and  Korpela  <  1983  )  both  performed  numerical 
studies  using  the  complete  Navier — Stokes  equations.  An  attempt 
was  made,  in  this  study,  to  duplicate  those  results  using  the 
boundary-layer  equations. 

Figure  27  shows  the  stream  function  contours  for  a  fluid  of 
Pr  —  1000.  de  Vahl  Davis  and  Mallinson  (  1975  )  reported 
multiple  ceil  formation  at  a  Rayleigh  number  between  2.4  x  105 
and  5.0  x  10J  .  They  based  their  Rayleigh  number  on  the  width 
of  the  vertical  slot  (  L  ).  The  Rayleigh  number  used  in  the 
boundary-i ayer  formulation  was  based  on  height  (  1  ).  An 
estimate  for  the  value  of  G  where  transition  should  occur  was 
made  by  converting  their  Rayleigh  number  in  the  following 
manner: 

Ra.  =  AT  (i  L3  — 

L  va 

Ra,  =  AT  /?  1 3  - — 

1  '  va 

G  =  L/l  . 


V 


10 


By  rearranging  the  last  equation  and  substituting  for  L, 


Ra, 


AT  ft  G313  2_ 
ua 


Rax  G' 


Now,  recalling  that 


~  1/4 

G  =  G  Raj  , 


gives  the  final  result  for  the  transformation, 

~  1/4 

G  =  (  G  RaL  )  . 

By  substituting  the  values  of  Ra^  =  2.4  x  10^  and  5.0  x  10^  ,  and 
G  =  0.1,  transition  was  expected  to  occur  between 

12.5  <  G  <  15.0  . 

Indeed,  multiple  cells  did  form  at  a  value  of  G  =  12.5  (  Figure 
27d  )  . 


Prior  to  transition,  the  stream  function  plots  agree  very 
well  with  the  results  obtained  by  de  Vahl  Davis  and  Mallinson 
and  Lee  and  Korpela  <  Figures  27a,  b,  and  c  ).  However,  once 
the  cells  formed,  the  arrangement  of  the  cell  patterns  generated 
by  the  boundary-1 ayer  program  <  Figures  27d  and  e  )  differed 
from  the  ones  that  they  obtained.  The  cell  patterns  that  they 
generated  showed  a  large  single  cell  in  the  center  of  the  box 
with  weaker  ones  just  above  and  below.  As  the  Rayleigh 
number  increased,  the  cells  became  stronger  and  weaker  cells 
began  to  form  between  them,  as  well  as  above  and  below. 
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The  boundary-layer  equations  generated  two  cells  of 
approx i matel y  the  same  strength  (  Figure  27d  )  ,  simultaneously. 

As  G  increased,  these  cells  became  stronger  and  began  to  split 
apart  to  form  other  cells.  Weaker  cells  also  formed  above  and 
bel ow. 

The  reason  for  the  differing  results  may  be  due  to  the 
boundary-1 ayer  approximation  where  the  streamwise  viscous  terms 
and  related  boundary  conditions  were  neglected.  For  finite 
Rayleigh  numbers  and  high  Prandtl  numbers,  these  terms  are  no 
longer  negligible  but  need  to  be  included.  However,  even 
without  these  terms,  multiple  cells  do  form;  but  in  a  different 
manner,  which  indicates  a  fairly  weak  dependency.  And,  in 
contrast  to  the  low  Prandtl  number  fluids,  the  fluid  motion  is 
now  slightly  skewed  in  the  opposite  direction.  This  also  agrees 
with  the  results  obtained  by  de  Vahl  Davis  and  Mallinson  as  well 
as  Lee  and  Korpela.  Finally,  it  should  also  be  noted  that  the 
solutions  obtained  for  the  high  Prandtl  number  fluids  are 
unsteady  in  nature. 

Figure  28  shows  the  vorticity  distribution  within  the  slot. 
Notice  that  for  high  Prandtl  number  fluids  the  vorticity  no 
longer  concentrates  in  distinct  nodal  regions  along  the  vertical 
walls  as  it  did  for  low  Prandtl  numbers;  but  instead, 
concentrates  almost  uniformly  along  the  entire  length.  In 
addition,  the  vorticity  distribution  becomes  anti  symmetrical 
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within  the  left  and  right  regions  (  Figures  28c,  d,  and  e  )  and 
no  longer  seems  to  be  influenced  by  the  number  of  convective 
cells  that  form  or  their  location  (  see  Figures  6,  12,  and  19  ). 

Temperature  distributions  (  Figure  29  )  obtained  from  the 
boundary  layer  program  also  agree  well  with  the  results  reported 
by  de  Vahl  Davis  and  Mallinson  and  Lee  and  Korpela,  up  to  the 
point  of  transition.  Notice  the  sinusoidal  appearance  of  the 
temperature  contours  prior  to  multiple  cell  formation  (  Figures 
29b  and  c  ) ,  in  contrast  to  the  nearly  vertical  contours  obtained 
for  the  low  Prandtl  number  fluids.  Notice  also  the  slight 
temperature  inversions  that  exist  where  cells  have  formed 
(  Figures  29d  and  e  ) .  Although  the  boundary-layer  program 
captured  these  temperature  inversions,  it  did  not  do  so 
completely.  It  failed  to  capture  the  temperature  inversion  at 
the  center  of  the  slot  because  no  cell  formed  in  that  region. 

And  finally,  analytical  (  Eqs  3.132  -  3.134  )  and  numerical 
comparisons  for  stream  function,  vorticity,  and  temperature  are 
shown  in  Figures  (  30-35  ).  Once  again,  the  comparisons  are 


very  good. 


Figure  3.  Stream  Function  contours  -for  Pr  =  0  . 

a)  G  =  3.0,  b)  6  =  3.6,  c>  G  =  3.7,  d>  G  =  4.0,  e)  G  =  4.2, 
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Figure  6.  Vorticity  Contours  -for  Pr 
a)  G  =  3.0,  b>  G  *  3.6f  c)  G  =  3.7 
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Figure  7.  Vorticity  variation  at  the  horizontal  centerline 
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Figure  9.  Stream  Function  contours  -for  Pr 
a)  S  =  1.0,  b)  6  =  1.5,  c)  S  *  1.6,  d) 
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centerline. 


PERCENT  ERROR 
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Figure  12.  Vorticity  Contours  -for  Pr 
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Figure  13.  Vorticity  variation  at  the  horizontal 
centerline. 
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Figure  18. 
a)  G  —  3.5, 


Stream  Function  contours  -for  Pr 
b)  G  =  3.9,  c)  G  =  4.0,  d) 
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Figure  21.  Stream  function  variation  at  the  horizontal 
centerline 
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Figure  23.  Vorticity  variation  at  the  horizontal  centerline. 
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Figure  25.  Temperature  variation  at  the  horizontal 
centerline 
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Figure  27. 
a )  G  =  3.0, 


Stream  Function  contours  for 
b)  G  =  10.0,  c>  G  =  12.0, 


Pr  =  1000  . 
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Figure  28.  Vorticity  Contour*  -for  Pr  *  10O0 
a)  G  =  3.0,  b)  6  «  10.0,  c)  G  =  12.0,  d> 


Figure  29.  Temperature  Contours  -for  Pr  =  1000  . 

a)  G  »  3.0,  b)  S  *  10.0,  c)  S  -  12.0,  d>  G  *  12.5,  e)  6  *  13.0 
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Figure  30.  Stream  function  variation  at  the  horizontal 
centerline. 
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Figure  31.  Stream  function  relative  error 
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Voriticity  variation  at  the  horizontal  centerline. 
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Figure  33.  Vorticity  relative  error. 


Figure  34.  Temperature  variation  at  the  horizontal 
centerline. 
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Figure  35.  Temperature  relative  error. 


VI.  Concl usi ons  and  Recommendations 


Concl usi ons. 

In  the  course  of  this  study,  it  was  found  that  there  are 
two  distinct  classes  of  fluid  behavior.  Fluids  with  low  Prandtl 
numbers  displayed  patterns  of  cell  formation,  vorticity,  and 
temperature  profiles  similar  to  each  other,  but  different  from 
those  with  high  Prandtl  numbers.  Below  Pr  =  0. 1  ,  thermal 
diffusion  dominates  and  the  instabilities  are  hydrodynamic  in 
nature.  Above  Pr  =  10.0  ,  momentum  diffusion  dominates  and  the 
instabilities  become  thermal  in  nature. 

The  2-D  boundary-layer  form  of  the  Navier — Stokes  equations 
were  able  to  capture  all  of  the  important  physics  of  multi¬ 
cellular  flow  in  the  vertical  slot  for  low  Prandtl  numbers  and 
Prandtl  number  of  order  one.  Excellent  qualitative  agreement 
with  results  obtained  by  other  researchers,  using  different 
methods,  was  achieved  with  this  approach.  Cells  formed,  gained 
in  strength,  and  distributed  themselves  throughout  the  vertical 
slot,  in  a  like  manner. 

For  very  low  Prandtl  numbers  (  ~  0.02  and  lower  ), 
convective  cells  formed  near  the  horizontal  boundaries  first, 
and  gained  in  strength  as  G  <  or  B,  as  the  case  may  be  ) 
increased.  The  weakest  cells  formed  in  the  center  of  the  slot. 
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The  opposite  was  true  for  air  (  Pr  =  0.71  ).  For  air,  the  cells 
formed  in  the  central  region  first,  and  were  strongest  there, 
while  weaker  cells  formed  near  the  horizontal  walls. 

For  low  Prandtl  numbers,  vorticity  was  found  to  have  a 
tendency  to  concentrate  in  nodal  regions  along  the  left  and 
right  vertical  boundaries.  Within  the  slot,  the  distribution  of 
vorticity  seemed  to  be  influenced  by  the  number  of  convective 
cells  that  formed  and  their  location.  The  direction  of 
vorticity  changed  from  counter — clockwise  at  the  left  and  right 
vertical  walls,  to  clockwise  in  the  center  of  the  slot. 

Temperature  profiles  obtained  for  air  (  Pr  ~  0.7  )  and 
mercury  (  Pr  ~  0.02  )  with  the  boundary-layer  approach  also 
showed  similar  trends  and  behavior  to  that  obtained  by  others. 
Distinct,  oscillatory  isotherms  were  present  in  the  multi¬ 
cellular  regions. 

In  all  cases,  analytical  results  were  in  excellent 
agreement  with  the  numerical  results  obtained  for  steady-state 
conditions. 

High  Prandtl  number  fluids  behaved  in  a  different  manner. 
The  boundary-layer  approach  was  not  able  to  completely  duplicate 
the  results  obtained  by  other  researchers.  This  is  most  likely 
because  the  streamwise  viscous  terms,  neglected  under  the 
assumptions  employed  in  the  boundary-1 ayer  technique,  have  a 
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more  pronounced  effect  on  cell  -formation  and  distribution  -for 
high  Prandtl  numbers  as  opposed  to  the  low  Prandtl  number  cases. 
None-the-1 ess ,  the  boundary-1 ayer  approach  was  able  to  capture 
single  and  multiple  convective  cell  formation  as  well  as  certain 
other  physical  properties  o-f  the  -flow.  Most  notable  among  these 
is  the  way  in  which  vorticity  and  temperature  develop  within  the 
slot . 


For  high  Prandtl  number  -fluids  (  Pr  ~  1000  )  ,  vorticity  no 
longer  concentrated  in  nodal  regions;  but  was  -fairly  uniform 
along  the  vertical  boundaries,  which  is  most  likely  a  ramifica¬ 
tion  of  large  momentum  di f f usi vi ties.  In  addition,  the 

vorticity  contours  appeared  to  be  antisymmetrically  distributed 
within  the  slot  and  no  longer  influenced  by  the  number  of  cells 
that  farmed  or  their  location. 

Temperature  profiles  were  highly  oscillatory  in  nature, 
even  before  multiple  cells  formed.  This  suggests  that  the 
instability  is  thermal  in  nature.  When  multiple  cells  did  form, 
the  temperature  profiles  showed  slight  inversions  present  in  the 
center  of  the  cells.  This  is  in  agreement  with  that  observed  by 
other  researchers.  However,  the  inversions  obtained  by  the 
boundary-1 ayer  method  were  not  as  pronounced. 

Overall,  the  boundary-1 ayer  approach  works  quite  well  for 
fluids  whose  Prandtl  number  is  of  order  one  and  lower  in  a 
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narrow  vertical  slot.  It  does  have  difficulty  modeling 
convective  cell  formation  at  high  Prandtl  numbers;  however,  the 
approach  captures  some  of  the  important  physical  properties. 
Although  quantitative  experimental  comparisons  are  lacking  in 
this  study,  comparisons  with  analytical  results,  as  well  as  the 
numerical  results  of  other  researchers,  lends  validity  to  this 
method. 

Recommendations. 

Additional  analytical  and  numerical  research  is  recommended 
with  the  high  Prandtl  number  fluids.  The  boundary-1 ayer  program 
can  be  modified  to  include  the  streamwise  viscous  terms 
neglected  in  this  study.  Incorporating  these  terms  would  also 
change  the  boundary  conditions  along  the  horizontal  walls.  The 
results  obtained  should  be  in  better  agreement  with  those 
obtained  in  previous  numerical  studies.  By  doing  this,  one  can 
demonstrate  the  importance  of  retaining  the  streamwise  viscous 
terms  for  large  Prandtl  number  fluids.  In  addition,  studies 
could  be  done  with  other  boundary  conditions  such  as  conducting 
horizontal  walls  or  horizontal  walls  with  constant  heat  flux. 

Due  to  insufficient  time,  the  Pr  -»  oo  case  was  not 
investigated.  This  study  could  be  done  also.  The  analytical 
and  numerical  expressions  for  the  limiting  case  of  Pr  -»  oo  are 
included  in  Chapters  III  and  IV. 
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Another  interesting  study  would  be  to  analyze  the  parallel, 
vertical  <  or  horizontal  )  -flat  plate  problem  using  the 
boundary-1 ayer  assumptions  employed  in  this  study. 

The  entropy  distribution  within  the  vertical  slot  could 
also  be  determined  and  platted.  This  could  be  accomplished  for 
the  entire  range  of  Prandtl  numbers  (  0  <  Pr  <  oo  )  and  would 
help  shed  more  light  on  the  multicellular  phenomenon. 

And  finally,  a  numerical  study  of  the  vertical  slot  with 
discrete  heat  sources  at  the  left  wall  could  be  performed.  As 
mentioned  in  the  literature  search  (  Chapter  I  ),  experimental 
data  exists  which  can  be  compared  to. 
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Appendix  A 


Supplementary  Analytical  Results 
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Figure  36.  Stream  function  variation  at  node  1 2. 
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Figure  38. 


Vorticity  variation  at  node  12. 
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Figure  39.  Vorticity  relative  error. 
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Figure  42.  Stream  function  variation  at  node  32. 
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Vorticity  variation  at  node  32. 
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Vorticity  relative  error. 
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Figure  46. 


Temperature  variation  at  node  32. 
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Figure  47.  Temperature  relative  error. 
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Figure  48.  Stream  function  variation  at  node  50. 
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Figure  49.  Stream  function  relative  error. 
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Vorticity  variation  at  node  50. 
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Vorticity  relative  error. 
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By  assuming  a  large  Rayleigh  number  and  small  gap,  the 
two-dimensional  Navier — Stokes  equations  are  reduced  to  the 
simpler  boundary-1 ayer  form.  Boundary-1 ayer  equations  are  then 
derived  for  both  the  low  (  Pr  i  0.1  )  and  high  <  Pr  >  10  ) 

Prandtl  number  regimes,  including  the  limiting  conditions  of 
Pr  -*  0  and  Pr  -»  oo.  Furthermore,  the  equations  obtained  for  the 
low  Prandtl  number  regime  also  hold  for  any  finite  Prandtl  number 
in  the  range  0  <  Pr  £  1.0. 

Numerical  solutions  to  these  equations  are  obtained  by 
using  a  stable,  second-order ,  f ul ly-impl ici t ,  time-accurate, 
Gauss-Si edel  iterative  procedure.  In  addition  to  numerical 
solutions,  steady-state  analytical  solutions  are  derived  using 
perturbative  methods. 

For  most  of  the  Prandtl  numbers  studied,  the  boundary-1 ayer 
approximations  produced  results  that  compared  well  with  those 

obtained  by  other  researchers  using  the  complete  Navier  Stokes 
equations  and  more  sophisticated  numerical  and  analytical 
techniques.  The  exception  is  in  the  high  Prandtl  number  regime, 
which  may  have  been  due  to  neglecting  the  streamwise  viscous 
terms.  However,  the  results  show  that  the  boundary-layer 
approach  captures  all  of  the  important  physics  of  low  Prandtl 
number  fluid  flows,  flows  where  Pr  ~  1.0,  and  a  significant 
portion  of  that  involved  with  fluids  of  high  Prandtl  number. 

Due  to  lack  of  time,  the  limiting  condition  where  Pr  ■*  oo 
was  not  accomplished.  However,  the  equations  for  this  type  of 
fluid  are  developed  and  presented  in  Chapter  III.  The  Pr  -*  » 
case  should  be  investigated  in  any  follow-on  studies  that  are 


conducted 


